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SECTION  1 


INTRODUCTION 


Design  and  evaluation  of  radio  frequency  systems  that  must  operate  through 
ionospheric  disturbances  resulting  from  high  altitude  nuclear  detonations  requires  an 
accurate  channel  model.  This  model  must  include  the  effects  of  scintillation  caused 
by  the  ionosphere  and  the  effects  of  high  gain  antennas  that  may  be  used  to  receive 
the  signals.  Such  a  model  can  then  be  used  to  construct  realizations  of  the  received 
signal  for  use  in  digital  simulations  of  trans- ionospheric  links  or  for  use  in  hardware 
channel  simulators. 

The  propagation  channel  is  conveniently  represented  in  terms  of  the  time- 
varying  channel  impulse  response  function  /i(r,f)  or,  equivalently,  by  its  Fourier  trans¬ 
form,  the  time-varying  channel  transfer  function  The  former  is  the  channel 

response  at  time  t  to  an  impulse  applied  at  t  -  T.  The  latter  is  the  channel  response 
to  a  sinusoidal  excitation  at  radian  frequency  w.  The  received  signal  u(t)  is  the  con¬ 
volution  of  the  channel  impulse  response  function  and  the  transmitted  modulation 
m(t): 


\ 

\ 


h(T,t)m(t  —  t)(Lt 


(i.i) 


The  purpose  of  this  report  is  to/describdHhe  Fortran  channel  model  ACIRF 
(Antenna/Channel  Impulse  Response  Function)  whiel)>generates  random  realizations 
of  the  impulse  response  function  at  the  outputs  of  multiple  antennas.  The  theory 
behind  the  impulse  response  function  has  been  discussed  by  numerous  authors  (see, 
for  example,  Knepp  1982  or  Dana  1986).  The  latter  reference  describes  the  theory 
of  propagation  through  random  media  and  antenna  aperture  effects  applications,  and 
is  one  of  the  references  from  which  ACIRF  was  developed.  The  other  reference  is 
Wittwer  and  Dana  (1989).  Channel  simulation  techniques  without  antennas  have  been 
described  by  Wittwer  (1980)  and  Knepp  and  Wittwer  ^i984).  The  former  reference 
contains  a  listing  of  the  original  CIRF  program  developed  by  Dr.  Leon  A.  Wittwer 
of  the  Defense  Nuclear  Agency.  Simulation  techniques  have  been  extended  to  include 
antenna  effects  in  Dana  (1986)  and  Wittwer  and  Dana  (1989).  ACIRF  is  intended 
replace  the  CIRF  code.  Thus  ACIRF  contains  the  channel  models  in  CIRF  without 
antennas  plus  it  has  new  channel  models  that  include  antenna  aperture  effects. 


Y- 
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Section  2  of  this  report  briefly  outlines  the  theory  behind  ACIRF.  The  Gen¬ 
eralized  Power  Spectra!  Density  (GPSD)  function,  which  describes  the  second  order 
statistics  of  the  channel  impulse  response  function,  is  discussed.  The  antenna  beam 
profile  parameters  (beamwidths  and  pointing  angles)  are  defined,  and  the  effects  of 
the  antenna  on  the  received  signal  statistics  are  described.  Channel  simulation  tech¬ 
niques  used  in  ACIRF  are  presented  in  Section  3.  Section  4  presents  examples  of  the 
output  of  a  matched  filter  which  illustrate  the  effects  of  scintillation  and  antennas  on 
the  received  signal.  The  inputs  and  outputs  of  the  ACIRF  Fortran  code  are  discussed 
in  Appendix  A. 
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SECTION  2 


THEORY 


The  starting  point  for  channel  modeling  is  the  generalized  power  spectral 
density  (GPSD).  This  function  describes  the  second-order  statistics  of  the  signal  in¬ 
cident  on  the  face  of  an  antenna.  When  scintillation  is  intense  or  fully  developed,  the 
first-order  signal  amplitude  statistics  are  well  described  by  the  Rayleigh  probability 
distribution,  and  the  signal  phase  is  uniformly  distributed  over  a  27t  radian  interval. 
The  strong  scattering  limit  can  also  be  obtained  by  applying  the  central  limit  theo¬ 
rem  to  the  superposition  of  the  many  randomly  scattered  waves  which  comprise  the 
received  signal  under  these  conditions.  The  in-phase  (/)  and  quadrature-phase  (Q) 
components  of  the  signal  are  then  independent,  zero-mean  Gaussian  random  variables 
with  variances  equal  to  one-half  the  total  signal  power.  These  conditions  are  both 
necessary  and  sufficient  for  Rayleigh  amplitude  statistics. 


2.1  GENERALIZED  POWER  SPECTRAL  DENSITY  (GPSD). 


The  derivation  of  the  GPSD  starts  with  Maxwell’s  equations  from  which 
the  parabolic  wave  equation  is  derived  (see,  for  example,  Dana  1986).  A  necessary 
condition  for  the  parabolic  wave  equation  to  be  valid  is  that  the  phase  perturbation 
over  a  distance  comparable  to  a  wavelength  be  small  compared  to  one  radian.  A 
sufficient  condition  is  that  the  angular  deviation  of  the  wave  relative  to  the  principal 
propagation  path  be  small  compared  to  one  radian.  These  conditions  are  generally 
satisfied  whenever  attenuation  of  the  propagating  wave  is  not  significant. 


2.1.1  Mutual  Coherence  Function. 

The  parabolic  wave  equation  can  be  solved  to  give  the  received  electric  field 
for  a  specific  distribution  of  the  index  of  refraction.  The  difficulty  is  that  the  index 
of  refraction  is  a  random  process.  The  parabolic  wave  equation  is  therefore  used  to 
derive  an  equation  for  the  two-position,  two-frequency,  two-time  mutual  coherence 
function,  T,  of  the  received  electric  field  E: 


r(z,y,u>,t)  =  (E(zi,yi,u>1,ti)E*(z2,y2,V2,t2))  (2.1) 
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where  x  =  x\  —  and  y  =  yi  -  jfe  are  relative  positions  in  the  plane  of  the  receiver, 
u  =  —  u>2  is  radian  frequency  difference,  and  t  =  ti  —  <2  is  time  difference. 

The  solution  of  the  differential  equation  for  the  mutual  coherence  function 
provides  a  description  of  the  second-order  statistics  of  the  received  signal.  The  time 
variation  part  of  the  mutual  coherence  function  takes  on  different  forms  depending 
on  the  nature  of  the  propagation  medium.  Two  limiting  models  are  of  interest.  The 
first  limit  is  the  “frozen-in”  model  wherein  the  striations  in  the  ionosphere  act  as  a 
cohesive  structure  that  moves  relative  to  the  propagation  path.  The  time  variations  of 
the  received  signal  are  then  a  straightforward  mapping  of  the  diffracted  electric  field 
as  it  drifts  past  the  receiver  antenna.  Thus,  in  this  frozen-in  model  there  is  a  strong 
coupling  between  the  spatial  and  temporal  variations  of  the  received  signal.  The 
second  limiting  model  is  a  “turbulent”  situation  wherein  striations  move  differently  at 
different  points  along  the  propagation  path,  causing  time  variations  of  the  diffracted 
electric  field  to  become  decorrelated  at  different  points  in  the  receiver  plane.  The 
turbulent  model  therefore  decouples  the  spatial  and  temporal  variations  of  the  received 
field. 


Wittwer  (1989)  has  recently  proposed  a  “general”  model  that  falls  between 
these  two  limiting  models.  This  model  is  currently  being  incorporated  into  a  new 
version  of  ACIRF  and  will  be  the  subject  of  a  future  volume  of  the  ACIRF  User’s 
Guide. 


The  mutual  coherence  function  for  the  frozen-in  model  is 


r(x,y,w,t) 
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The  coupling  between  spatial  and  temporal  variations  for  the  frozen-in  model  is  a 
result  of  the  x  —  £xf /r0  term.  For  the  turbulent  model,  the  mutual  coherence  function 
is 


T{x,y,u,t)  =  rT(t) 
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where  IV  (f)  is  the  temporal  coherence  function  which  does  not  depend  on  either  the 
spatial  variables  or  frequency. 

By  examining  these  two  forms  for  the  mutual  coherence  function,  it  can 
be  seen  that  the  decorrelation  distance  ts  is  the  1/e  point  of  r(x, 0,0,0),  and  the 
decorrelation  distance  ly  is  the  1/e  point  of  r(0,y,0,0).  The  asymmetry  factor  A  is  a 
function  of  the  two  decorrelation  distances: 


U  +  (4/4)4J 

The  coordinate  system  of  the  scattering  is  defined  by  the  line  of  sight  di¬ 
rection  z  and  the  geomagnetic  field  B  at  the  altitude  in  the  ionosphere  where  the 
scattering  occurs.  The  orientation  of  the  x  axis, 

.  _  B  x  z 
X  ~  |JB|  sin  $  ’ 

is  chosen  such  that  lx  <  £„.  This  scattering  coordinate  system  is  shown  in  Figure  2.1. 
The  penetration  angle  $  is  the  angle  between  the  line  of  sight  direction  and  the 
magnetic  field. 

In  general,  both  decorrelation  distances  should  be  part  of  the  specification 
of  the  channel  parameters.  However,  if  only  the  minimum  decorrelation  distance  ix 
(also  referred  to  as  £o  in  the  literature)  is  specified,  then  the  ratio  of  the  decorrelation 
distances  £„/£*  can  be  estimated  from  the  formula: 


j-  =  ’\J cos2  $  +  <72  sin2  $  (2.5) 

where  q  is  the  axial  ratio  which  is  equal  to  the  ratio  of  the  characteristic  length  of 
striations  along  the  field  lines  to  the  characteristic  thickness  in  a  direction  normal  to 
the  field  lines.  The  standard  value  for  q  has  been  15.  However,  recent  research  on  the 
physics  of  striations  indicates  that  q  may  be  as  large  as  70  (Wittwer  1988). 

The  decorrelation  time  r0  is  the  1/e  point  of  r(0,0,0,  t)  for  the  frozen-in 
model.  In  order  to  maintain  a  consistent  definition  of  decorrelation  time,  it  will  also 
be  required  that  rr(r0)  equal  1/e  for  the  turbulent  model. 
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Line  of  Sight  Penetration 


Direction  of  Minimum 
Decorrelation  Distance 

Figure  2.1.  Scattering  coordinate  system. 

The  coherence  bandwidth,  wCOfc,  is  defined  using  the  GPSD  which  will  be 
discussed  next.  The  quantity  a  is  a  delay  parameter  which  is  equal  to  the  ratio  of 
the  time  delay  spread  caused  by  diffraction  to  that  caused  by  refraction.  In  strong 
scattering  where  diffractive  effects  dominate,  the  value  of  a  is  fairly  large. 


2.1.2  Angular-Delay  Power  Spectral  Density. 

The  Fourier  transform  of  the  mutual  coherence  function  is  the  generalized 
power  spectral  density  of  the  received  signal: 

/DO  fOO  f  00  (hj  [  00 

dx  iy  /  —  dtT(x,y,u,t)  x 

-oo  J  -oo  J  —  oo  Z7C  J~oo 

exp  [-i(Kxx  +  Kvy  -  u/r  +  wpf)]  (2.6) 
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where  up  is  the  Doppler  radian  frequency,  r  is  the  time  delay,  and  Kx  and  Ky  are 
related  to  the  scattering  angles  6X  and  0V  about  the  x  and  y  axes  as  follows: 


2n  sin  6X 

K-  =  ~r~ 

^  _  2n  sin  6V 

Kv  -  j 

The  quantity  S(KX,  Kv,  T,up)dKxdKydrdup  is  the  mean  signal  energy  incident  on  the 
face  of  a  receiving  antenna  arriving  from  angles  Kx  to  Kx  +  dKx  and  Kv  to  Ky  +  dKv\ 
with  delays  in  the  interval  r  to  r  +  dr,  and  with  Doppler  frequencies  in  the  interval 
Up  to  Up  +  duip. 

In  general,  the  GPSD  can  be  written  as  a  product  of  a  Doppler  term  times 
an  angular-delay  term, 


(2.7) 


S(K„K„r,uB)  =  SD(wD)SKr(K„K„,r)  , 


(2.8) 


where  the  angular-delay  part  of  the  GPSD  is 


Sk,(Kx,K„t)  = 


7T 

L2J 


1/2 


aoJcohZzty  exp 


KX  +  KX 


exp 


a" 


^coh^ 


K(Kl  +  KX 


l2' 


(2.9) 


2.1.3  Doppler  Spectra. 

The  Doppler  part  of  the  GPSD  has  different  forms  for  the  frozen-in  and 
turbulent  models.  For  the  frozen-in  model,  the  drift  velocity  of  striations  is  chosen  to 
be  along  the  x-direction  in  the  receiver  plane.  The  Doppler  component  of  the  GPSD 
is  then 


Sp(up)  =  2n6  (upTo  —  Kxlx) 


(2.10) 
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where  £(•)  is  the  Dirac  delta  function.  This  choice  for  the  orientation  of  the  drift 
velocity  minimizes  the  value  of  To.  The  decorrelation  time  in  the  frozen-in  model  is 
the  time  required  for  the  diffraction  pattern  to  drift  one  decorrelation  distance  in  the 
receiver  plane. 

Combining  Equations  2.8  through  2.10  and  integrating  over  angle  and  delay 
gives  the  Doppler  spectrum  for  the  frozen-in  model: 


S(wD) 


°°  dKx  [«>  dKy 


r°°  mu  r 

J-oo  2n  J- 


2tt 


\pn  tq  exp  — 


2  2  " 
TqU}) 


f  dr  Sd{^d)Skt{Kx,Kv,t) 

J -oo 

.  (Frozen-In  Model) 


(2.11) 


The  temporal  coherence  function  for  this  case  also  has  a  Gaussian  form: 


rT(t)  =  f  exp  {iuDt)S(uD)^^  =  exp 

J —oo  &7T 


(2.12) 


For  the  turbulent  mode’ ,  the  Doppler  component  of  the  GPSD  is  a  function 
only  of  the  Doppler  frequency.  In  this  model  there  is  no  coupling  between  scattering 
angle  and  Doppler  frequency  or,  equivalently,  between  position  and  time.  The  Doppler 
spectrum  is  therefore  independent  of  delay  so  the  temporal  variations  exhibit  the  same 
rate  at  all  delays.  A  mathematically  convenient  form  of  the  Doppler  spectrum  for  the 
turbulent  model  is  the  so-called  f~A  spectrum: 


S(u>d)  = 


4  Tq _ 1 _ 

°  [l  +  7o w|)/o2]2 


(Turbulent  Model) 


(2.13) 


where  the  parameter  a  (a  =  2.146193)  is  determined  by  the  condition  that  IV(ro) 
equal  1/e.  The  corresponding  temporal  coherence  function  then  has  the  form: 


2.1.4  Frequency  Selective  Bandwidth  and  Angle  of 
Arrival  Variance. 

An  important  measure  of  the  effects  of  scintillation  on  the  propagation  of 
wide  bandwidth  signals  is  the  frequency  selective  bandwidth.  This  quantity  is  defined 
as 


/o  = 


1 

27 rcy 


(2.15) 


where  oT  is  the  standard  deviation  of  the  time  delay  jitter  of  the  received  signal: 


o*  =  (r2)  ~  {r? 


(2.16) 


The  moments  of  delay  are  calculated  using  the  angular-delay  part  of  the  GPSD, 


<r">=r  srsr  ***  *»«.*..'> 


r°°  dKu  r°° 


(2.17) 


These  equations  then  give  the  relationship  between  ueoh  and  the  frequency  selective 
bandwidth  /0: 


^eoh  —  2tt/o\/i  +  1/a2  .  (2.18) 

An  important  signal  parameter  for  systems  that  utilize  large  aperture  re¬ 
ceiving  antennas  is  the  angular  scattering  variance.  Whenever  the  striations  in  the 
ionosphere  are  not  aligned  with  the  propagation  direction,  which  is  typically  the  case, 
the  angular  scattering  will  be  different  about  the  x  and  y  axes.  The  standard  devia¬ 
tions  of  the  angle  of  arrival  fluctuations  about  the  x  and  y  axes  are: 


Oix  = 


Ogy  — 


A 

Vlvlz 

A 

y/2n£v 


(2.19) 
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If  the  propagation  path  is  aligned  with  the  striations  in  the  ionosphere,  the 
decorrelation  distances  in  the  x  and  y  directions  will  be  the  same.  In  this  isotropic 
case  the  scattering  angles  are  the  same  in  each  direction.  A  three-dimensional  plot 
of  the  isotropic  one-dimensional  GPSD  is  shown  in  Figure  2.2.  This  plot  shows  an 
example  of  the  mean  received  power  as  a  function  of  normalized  angle-of-arrival,  KIq, 
and  normalized  time  delay,  u^r.  For  this  isotropic  case,  the  quantity  K  is  either  Kx 
or  Kv  since  they  are  equal.  Similarly,  the  quantity  Iq  is  either  tx  or  iv.  The  value  of  a 
for  this  figure  is  taken  to  be  4,  so  consequently  the  quantity  ucoti  is  essentially  equal  to 
27r/0.  It  can  be  seen  from  the  figure  that,  the  power  arriving  at  large  angles  is  also  the 
power  arriving  at  long  delays.  The  power  arriving  at  long  delays  thus  is  more  likely 
to  be  attenuated  if  the  receiving  antenna  has  a  narrow  beamwidth. 


Figure  2.2.  Example  of  the  generalized  power  spectral  density. 
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2.1.5  Diffraction  Limited  Form  of  GPSD. 


Under  the  strong  scattering  conditions  in  the  ionosphere  that  cause  signal 
scintillation  at  radio  frequencies,  diffractive  effects  dominate  refractive  effects.  In  this 
case,  the  value  of  a  in  the  GPSD  is  large.  In  the  limit  that  a  approaches  infinity,  the 
form  of  the  GPSD  can  be  simplified  by  noting  that 

cH^oo  [aexp(-a2x2)]  =  .  (2.20) 

In  this  limit,  the  angular-delay  part  of  the  GPSD  becomes 


Sgr  (F^z  ,  Ky  ,  7")  TTWjo/jfjfy  exp 


Kjll  +  Kl(l 


^coh^ 


a(jt|  +  Kl)il 


(2.21) 


This  geometric  optics  limit  results  in  a  delta  function  relationship  between  angle  and 
delay.  This  form  of  the  angular-delay  part  of  the  GPSD  will  be  used  to  generate 
realizations  of  the  impulse  response  function  at  the  outputs  of  antennas. 


2.2  ANTENNA  EFFECTS. 


The  GPSD  of  the  signal  at  the  output  of  an  antenna,  S^(K,t,ud),  is  given 
by  (Frasier  1988) 


Sa{K,t,ud)  =  G(K  -  K0)S{K,t,ud)  (2.22) 

where  G(K  —  Kq)  is  the  two  dimensional  beam  profile  of  an  antenna  pointing  in  the 
direction  K0.  The  effect  of  the  antenna  then  is  to  attenuate  the  received  signal  energy 
at  angles  outside  of  the  main  beam  of  the  antenna.  If  the  antenna  is  pointed  along 
the  line  of  sight,  the  signal  energy  that  is  attenuated  by  the  antenna  is  that  with  large 
angles  of  arrival.  This  is  also  the  energy  with  large  times  of  arrival.  The  effect  of  the 
antenna  will  then  be  to  reduce  the  received  signal  power,  and  also  to  reduce  the  delay 
spread  of  that  energy. 
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2.2.1  Antenna  Coordinate  System. 


The  antenna  coordinate  system  is  shown  in  Figure  2.3.  It  is  assumed  that 
the  face  of  the  antenna  is  in  the  x—y  plane  and  the  boresight  is  pointed  along  the  line  of 
sight.  The  beam  in  this  figure  is  pointed  away  from  the  line  of  sight  in  the  direction  K0 
defined  by  a  pointing  angle  tp  and  an  azimuth  angle  e.  The  rotation  angle  is  the  angle 
between  the  scattering  x  axis  and  an  antenna  axis  u.  The  antenna  (u,  v)  coordinate 
system,  where  v  is  in  the  x  —  y  plane  and  is  orthogonal  to  u,  is  chosen  for  convenience 
in  describing  the  antenna  beam  profile.  For  example,  if  the  antenna  is  rectangular, 
then  the  (u,  t))  axes  should  be  aligned  with  the  sides  of  the  aperture.  In  this  coordinate 
system,  the  antenna  beam  profile  G(0U,0V )  is  approximated  as  a  product  G(0U)  times 
G(0V )  where  0U  and  0V  are  angles  about  the  u  and  t;  axes  respectively.  In  many  cases, 
the  axis  u  will  be  parallel  to  the  local  earth  tangent  plane. 

The  components  of  the  antenna  pointing  direction  K0  in  terms  of  the  rota¬ 
tion,  pointing,  and  azimuth  angles  are 


2tt 

K0x  =  —  sin  rp  cos(x  +  e) 

(2.23) 

K0v  =  —  sin  rp  sin(x  +  e)  . 

A 

(2.24) 

2.2.2  Antenna  Beam  Profile. 

The  antenna  beam  profile  is  assumed  to  have  a  Gaussian  shape  and  to  be 
separable  in  the  u  —  v  coordinate  system: 

G(K,  K.)  =  exp  \-a\Kl  -  .  (2.25) 

The  peak  gain  G(0,0)  has  been  set  to  unity  because  this  value  is  usually  included  in 
the  calculation  of  the  mean  received  power.  For  either  the  u  or  the  v  direction,  the 
antenna  beam  profile  can  also  be  written  as 


G{0)  =  exp[-(ln2)(20/0o)2]  (2.26) 

where  60  is  the  3  dB  beamwidth  [i.e.  full  width  at  half  maximum,  G(0o/2)  =  1/2). 
Equating  these  two  profiles  gives 
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Line  of  Sight 

z 


(2.27) 


a 


2 

£ 


(ln2)X2 


where  0O$  is  the  3  dB  beamwidth  in  the  £  =  u  or  the  £  =  v  direction. 

In  the  next  subsections,  the  3  dB  beamwidths  are  related  to  antenna  sizes 
for  the  special  cases  of  uniformly  weighted  circular  and  rectangular  apertures. 


Uniformly  Weighted  Circular  Apertures 

For  a  uniformly  weighted  circular  aperture,  the  well  known  form  for  the 
power  beam  profile  is 


G(»)  = 


4J2[n[D/X)  sin^] 
[n{D/X)  sin0]J 


(2.28) 


where  Ji(-)  is  the  Bessel  function.  The  3  dB  full  width  at  half  maximum  beamwidth, 
in  terms  of  the  diameter  D ,  is 


0o  =  1.02899(A/D)  radians  .  (2.29) 

If  this  beam  profile  is  approximated  by  a  Gaussian  profile  with  the  same  3  dB 
beamwidth,  then  the  coefficients  that  appear  in  Equation  2.24  are 


2  2  J/n2jI>2_ 

u  “  (1.02899 7r)2 

Uniformly  Weighted  Rectangular  Apertures 

For  uniformly  weighted  rectangular  apertures,  the  antenna  beam  profile  in 
either  the  £  =  u  or  £  =  v  direction  has  the  familiar  sinJ(x)/i2  form: 


sin2[7r(ZVA)sinflf| 
[tt(ZV  A)  sin  0{]2 


(2.31) 


where  D {  is  the  length  of  the  aperture.  The  3  dB  beamwidth  is 
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0O(  =  0.885893(A/I>4)  radians  . 


(2.32) 


and  the  coefficients  that  appear  in  Equation  2.24  are 

_  (ln2)D\ 

*  (0.885893tt)1  ’ 


(2.33) 


2.2.3  Antenna  Filtering  Equations. 

Frasier  (1988)  has  derived  filtering  equations  for  an  antenna  pointing  in  an 
arbitrary  direction  away  from  the  line  of  sight.  The  equations  for  the  scattering  loss, 
frequency  selective  bandwidth,  decorrelation  time  and  distances  are  reproduced  here. 
Further  details  on  the  derivation  of  these  equations  may  be  found  in  Frasier  (1988). 

For  an  arbitrary  antenna  with  a  beam  profile  characterized  by  two  beamwidths 
and  for  general  anisotropic  scattering,  it  is  mathematically  convenient  to  write  the 
angular-delay  part  of  the  GPSD  using  matrix  notation: 


SKt{K ,  r)  =  7 Tucohltly  exp\-K  •  L  •  K}6 


WcohT 


A  KHl' 
4 


(2.34) 


where 


K  =  jK[+lq  (2.35) 

is  the  magnitude  of  the  vector  K ,  and  the  decorrelation  length  matrix  L  in  the  x  —  y 
coordinate  system  is 


The  mean  received  power  at  the  output  of  the  antenna,  is 


(2.36) 
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(2.37) 


Pa  = 


det(L) 


det(Q) 


1/1 


G(K0)  exp(m  •  Q~l  •  m) 


where  G(K0)  is  the  line  of  sight  gain  of  the  antenna  pointed  in  the  direction  K0  away 
from  the  line  of  sight, 


Q  —  L  *1*  ct 


(2.38) 


m  —  a  -  K0  ,  (2.39) 

and  where  the  matrix  a,  in  the  x  —  y  coordinate  system,  is 

a2cos2x  +  a2sin2x  (a2  -  a2)  sin  x  cos  x 

a  = 

(a2  -  a2)  sin  x  cos  x  <*l  sin2  X  +  aj  cos2  x 

Note  that  Equation  2.36  is  independent  of  the  choice  of  coordinate  systems  as  long  as 
all  of  the  matrices  are  expressed  in  the  same  coordinate  system.  The  scattering  loss, 
in  decibels,  of  the  antenna  is  just  -10  logi0(P/t)* 

The  antenna  filtered  values  of  the  decorrelation  distances  are 


(2.40) 


Iaz  =  2 


det{Qyl/2 


Q 


vv 


t-Ay  —  2 


det{Qy1/2 


Q 


xz  J 


(2.41) 


(2.42) 


where  Qxx  and  Qvv  are  the  two  diagonal  components  of  Q  in  the  x  —  y  coordinate 
system.  Note  that  the  filtered  decorrelation  distances  depend  only  on  the  beamwidths 
of  the  antenna  and  are  independent  of  pointing  angle.  However,  as  a  beam  is  pointed 
away  from  the  line  of  sight,  the  beamwidth  will  increase  because  of  foreshortening  of 
the  aperture.  Thus  the  independence  of  the  antenna  filtered  decorrelation  distances 
with  pointing  angle  is  only  true  for  small  pointing  angles. 
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The  antenna  filtered  decorrelation  time  for  the  frozen-in  model,  when  the 
relative  motion  of  the  antenna  and  random  diffraction  pattern  is  in  the  x  direction,  is 


Ta  _  Ux 
Tq  tx 


2~ 

in 


det(Q) 


Q 


vv  J 


1/2 


(Frozen-in  Model) 


(2.43) 


The  expression  for  ta/tq  is  also  independent  of  pointing  angle.  For  the  turbulent 
model, 


—  =  1  .  (Turbulent  Model)  (2-44) 

TO 


Finally,  the  antenna  filtered  frequency  selective  bandwidth  is 


/a  _  1^ 

fo  det(L) 

The  form  of  this  result  emphasizes  its  invariance  with  respect  to  the  choice  of  coordi¬ 
nate  systems. 


tr(L') 


4m  •  Q~z  •  m  + 

det2(Q) 


(2.45) 


2.2.4  Antenna  Filtering  Effects. 

The  effect  of  antenna  filtering  is  illustrated  in  Figure  2.4,  which  shows  four 
plots  of  the  GPSD  at  the  outputs  of  four  uniformly  weighted  circular  antennas  for 
isotropic  scattering  (lx  =  ty  =  to) .  The  antennas  in  this  example  are  all  pointed  along 
the  line  of  sight.  The  upper  left  plot  is  the  same  as  in  Figure  2.2  and  is  for  a  point 
antenna  (D  to) ,  for  which  there  is  no  antenna  filtering.  The  other  three  plots  are 
for  cases  where  the  ratio  of  antenna  diameter  to  decorrelation  distance  is  equal  to  1, 
2,  and  4.  As  the  antenna  size  increases  for  a  given  value  of  to  or,  equivalently,  as  the 
decorrelation  distance  decreases  for  a  given  antenna  diameter  D,  more  of  the  energy 
arriving  at  large  scattering  angles  is  filtered  out. 

Antenna  filtering  has  two  major  effects,  one  of  which  is  generally  bad  and  the 
other  potentially  good  in  terms  of  system  performance.  The  bad  effect  is  scattering 
loss.  Energy  that  arrives  at  large  angles  is  not  captured  by  the  relatively  narrow  beam 
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Generalized  power  spectra  at  the  outputs  of  circular 
antennas. 


of  the  large  aperture  antenna.  An  equivalent  statement  is  that  spatial  decorrelation 
of  the  received  electric  field  across  the  antenna  aperture  causes  a  reduction  in  the 
effective  gain  of  the  antenna.  The  magnitude  of  the  scattering  loss  is  a  function  of 
the  ratio  of  the  antenna  diameter  to  decorrelation  distance  of  the  incident  wave  and 
a  function  of  the  pointing  angle  of  the  antenna  beam  away  from  the  line  of  sight. 

Figure  2.5  shows  the  scattering  loss  of  a  uniformly  weighted  circular  antenna 
versus  the  ratio  to/D  when  the  pointing  angle  of  the  antenna  beam  is  0,  60/2  and  60. 
It  is  assumed  for  this  example  that  the  beamwidth  remains  constant  as  the  beam  is 
pointed  away  from  the  line  of  sight.  When  the  decorrelation  distance  is  large  compared 
to  the  antenna  diameter,  the  curves  approach  the  line  of  sight  loss  of  beam  (0  dB  for  0 
pointing  angle,  3  dB  for  60/2  pointing  angle,  and  12  dB  for  $0  pointing  angle).  When 
the  decorrelation  distance  is  small  compared  to  the  antenna  diameter,  the  angle  of 
arrival  spread  of  the  incident  energy  is  large  compared  to  the  beamwidth,  so  the 
scattering  loss  is  insensitive  to  pointing  angle  as  long  as  the  pointing  angle  remains 
within  a  few  beamwidths  of  the  line  of  sight. 


Figure  2.5.  Scattering  loss  for  a  uniformly  weighted  circular  antenna. 
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The  potentially  beneficial  effect  of  antenna  filtering  is  that  the  energy  which 
is  filtered  out  is  that  which  arrives  at  relatively  large  delays.  It  is  this  delayed  energy 
that  causes  most  of  the  intersymbol  interference  in  the  detection  and  demodulation  of 
wide  bandwidth  signals.  Stated  another  way,  because  the  antenna  filters  out  much  of 
the  delayed  energy,  and  since  the  frequency  selective  bandwidth  is  an  inverse  measure 
of  the  signal  delay  spread,  the  filtering  increases  the  frequency  selective  bandwidth 
of  the  output  signal  relative  that  at  the  antenna  input.  This  is  shown  in  Figure  2.6, 
which  has  plots  the  ratio  of  the  filtered  to  unfiltered  frequency  selective  bandwidth, 
/a/ /o>  as  a  function  of  the  ratio  Iq/D  for  a  uniformly  weighted  circular  antenna. 

Figure  2.6  illustrates  an  interesting  effect  of  pointing  the  beam  away  from 
the  line  of  sight.  As  the  pointing  angle  increases  at  a  fixed  value  of  the  ratio  Iq/D, 
the  standard  deviation  of  the  time  delay  jitter  increases  and  the  filtered  value  of  the 
frequency  selective  bandwidth  decreases.  Thus  intersymbol  interference  effects  in  the 
output  of  a  beam  may  increase  as  the  beam  is  pointed  away  from  the  line  of  sight. 

Antenna  filtering  can  also  affect  the  decorrelation  time  of  the  received  signal 
as  shown  in  Figure  2.7.  The  filtered  value  of  decorrelation  time  is  independent  of 
pointing  angle,  but  depends  strongly  on  the  model  for  the  temporal  fluctuations.  For 
turbulent  model,  the  Doppler  spectrum  is  independent  of  angle,  so  the  signal  energy 
arriving  at  all  angles  has  the  same  decorrelation  time.  Hence  the  filtered  value  of 
decorrelation  time  is  equal  to  that  of  the  incident  signal  (r^/r0  =  l)  for  this  model. 
For  the  frozen-in  model,  there  is  strong  coupling  between  the  angular  spectrum  and 
the  Doppler  spectrum  of  the  incident  signal  energy.  The  effect  of  an  antenna  is  to 
narrow  the  angular  spectrum  the  received  signal  energy  as  shown  in  Figure  2.4.  Thus, 
for  the  frozen-in  model,  the  antenna  also  narrows  the  Doppler  spectrum  of  the  received 
signal  energy.  Because  the  decorrelation  time  of  the  signal  is  an  inverse  measure  of 
the  width  of  the  Doppler  spectrum,  the  decorrelation  time  of  the  signal  out  of  the 
antenna  increases  as  antenna  filtering  effects  increase.  Reality  should  lie  somewhere 
between  these  two  limiting  models. 
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Figure  2.6. 


Figure  2.7. 


Filtered  frequency  selective  bandwidth  for  a  uniformly 
weighted  circular  antenna. 


Filtered  decorrelation  time  for  a  uniformly  weighted  circular 
antenna. 
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SECTION  3 


CHANNEL  SIMULATION 


The  purpose  of  this  section  is  to  describe  statistical  channel  simulation  tech¬ 
niques  that  allow  realizations  of  the  impulse  response  functions  to  be  generated  at  the 
outputs  of  multiple  antennas  with  spatial  and  temporal  correlation  properties  given 
by  the  GPSD  and  with  Rayleigh  amplitude  statistics.  The  realizations  generated  with 
these  techniques  represent  only  the  diffractive  part  of  the  received  voltage  and  are 
valid  under  strong  scattering  conditions  where  the  GPSD  is  valid  and  where  Rayleigh 
statistics  apply.  Under  these  conditions,  the  impulse  response  function  represents  a 
solution  of  Maxwell’s  equations  for  propagation  of  RF  waves  through  randomly  struc¬ 
tured  ionization. 

The  basic  formalism  to  generate  statistical  realizations  of  the  channel  im¬ 
pulse  response  function  without  antenna  effects  explicitly  included  was  developed  by 
Wittwer  (1980)  for  isotropic  irregularities  and  was  extended  by  Knepp  (1982)  to  the 
case  of  elongated  irregularities.  (The  elongated  case  corresponds  to  a  90°  penetration 
angle  and  to  an  infinite  axial  ratio.)  Channel  simulation  techniques  have  been  gener¬ 
alized  in  ACIRF  to  include  the  effects  of  general  anisotropic  scattering  and  multiple 
receiving  antennas. 

Another  extension  of  previous  channel  modeling  techniques  that  is  used  in 
ACIRF  is  the  use  of  the  diffraction  limited  form  of  the  GPSD  (Equation  2.21).  In 
this  limit,  there  is  a  delta  function  relationship  between  scattering  angle  and  delay. 
This  allows  the  angular-delay  part  of  the  random  signal  to  be  generated  in  angle  only. 
The  relationship  between  angle  and  delay  is  then  used  to  associate  angular  spectral 
components  within  an  angular  annulus  to  a  unique  delay  bin.  The  advantage  of  this 
procedure  is  that  it  requires  one  less  dimension  in  the  arrays  used  for  channel  modeling 
thereby  reducing  both  computation  time  and  executable  code  size. 

Channel  simulation  techniques  rely  on  the  assumptions  that  the  channel 
is  statistically  stationary  and  that  the  first  order  statistics  of  the  impulse  response 
function  are  normal.  The  first  step  in  channel  modeling  is  to  generate  the  random 
angular  spectral  components  of  the  channel  before  the  antennas.  Because  the  channel 
is  stationary  in  space,  these  components  are  uncorrelated  from  one  angle  to  another. 
The  angular  spectrum  of  the  channel  is  multiplied  by  the  antenna  beam  profiles  to 
produce  the  angular  spectrum  at  the  outputs  of  the  antennas.  The  delta  function 
relationship  between  angle  and  delay  is  used  to  associate  the  angular  components  with 
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unique  delay  bins.  After  this  point,  techniques  used  for  the  frozen-in  and  turbulent 
models  to  create  signal  temporal  variations  are  quite  different. 

For  the  frozen-in  model,  temporal  variations  result  from  a  drift  of  a  ran¬ 
dom  diffraction  pattern  across  the  receiving  antennas.  Of  course  it  is  equally  valid 
to  assume  that  the  antennas  drift  past  the  random  diffraction  pattern.  The  latter 
viewpoint  is  used  in  the  channel  model.  The  impulse  response  functions  are  generated 
by  fast  Fourier  transforming  the  angular  spectra  at  the  outputs  of  the  antennas  to 
a  sequence  of  positions  along  the  x  axis.  Temporal  variations  result  by  setting  the 
relative  speed  between  the  diffraction  pattern  and  the  antennas  to  one  decorrelation 
distance  per  decorrelation  time  (i.e.  speed  is  equal  to  lx/r0).  There  is  then  a  one-to- 
one  relationship  between  the  Kx  angular  spectrum,  modified  by  Kx  component  the 
antenna  beam  profile,  and  the  Doppler  spectrum  as  described  by  Equation  2.9.  The 
output  of  the  channel  model  is  the  impulse  response  functions  h(jAr,nAt)  for  each 
antenna  at  a  set  of  discrete  times  nAf  and  delays  j'Ar. 

For  the  turbulent  model,  temporal  variations  in  the  impulse  response  func¬ 
tion  result  from  temporal  variations  in  the  random  angular  spectral  components.  An 
f~A  Doppler  spectrum  has  been  chosen  for  computational  convenience.  It  is  then  pos¬ 
sible  to  generate  the  temporal  variations  of  the  angular  spectral  components  directly 
in  the  time  domain  by  passing  white  Gaussian  noise  through  two-pole  low  pass  filters. 
At  a  given  time,  the  angular  spectral  components  of  the  channel  are  modified  by  the 
antenna  beam  profiles  and  divided  into  delay  bins.  Then  discrete  Fourier  transforms 
are  performed  from  the  angular  domain  to  the  locations  of  the  antennas.  At  the  next 
time  step,  the  angular  spectral  components  of  the  channel  are  updated  using  the  /~ 4 
filters,  modified  by  the  antenna  beam  profiles,  divided  into  delay  bins,  and  Fourier 
transformed  to  the  antenna  positions,  and  so  on.  Again  the  output  is  the  impulse 
response  functions  for  each  antenna  at  a  set  of  discrete  times  and  delays. 


3.1  FROZEN-IN  MODEL. 

The  steps  used  to  generate  realizations  of  the  impulse  response  functions 
at  the  outputs  of  multiple  antennas  using  the  frozen-in  model  are  outlined  below. 
Further  details  of  the  algorithm  may  be  found  in  Wittwer  and  Dana  (1989). 
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3.1.1  Discrete  Evaluation  of  the  GPSD. 


The  first  step  in  generating  the  impulse  response  function  is  the  evaluation 
of  the  GPSD  on  a  discrete  angular-delay  ( Kx  —  Ky  —  r)  grid.  The  angular-delay  grid 
centers  are  defined  by 


Kx  =  kAKz 
Ky  =  lAKy 

t  =  (j  +  1/2)A  r  (3.1) 

where  k ,  /,  and  j  are  integers.  The  angular  and  delay  grid  sizes  and  the  requirements 
on  the  number  of  grid  points  will  be  discussed  later.  For  the  jth  delay  bin,  the  GPSD 
will  be  non-zero  only  for  the  angles  that  satisfy  the  condition 

jucohAr  <  — * - - - - —  <  (j  +  l)ucoflAr  (3.2) 

4 

where,  for  this  diffraction  limited  case,  <*>„„>,  is  equal  to  2nfo- 

The  mean  signal  energy  in  a  Kx  —  Kv  grid  cell  at  the  output  of  the  mth 
antenna  is 


E„ 


r(k+l/2)bKz  dK  /■(r+l/3)AA’-  dK 

(M)=/  ~~  ^  Gm(Kx  -  K0x,Kv  -  K0v)Sk(Kx,Kv) 

'  Jk-1/2)&K,  2n  J(l-l/2)AK,  27T 


(3.3) 


where  the  angular  power  spectral  density  Sr  is 


Sk(Kx,Kv)  =  f  SKt{Kx,Kv,t)(It  =  ntxly  exp 
J -00 


KX  +  Kffl 


(3.4) 


Equation  3.3  is  completely  general,  but  it  implies  that  Em(k,  l)  must  be  computed  and 
stored  for  each  separate  antenna.  If,  however,  it  is  assumed  that  the  antenna  beam 
profile  is  constant  across  a  Kx  -  Ky  bin,  then  Equation  3.3  can  be  approximated  by 
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EJk,l)  =  Gm(kAKx  -  K„„IAK,  -  K„„)E(k,l) 


(3.5) 


where 


(*+l/2)A Kz  dK,  r(l+l/2)*K,  dKu 


E(k,l)  =  [  ’  [  ' 

J(k-1/2)AKI  2%  J(l- 


(Jt-l/2)Aff,  27T  J(l-l/2)AK,  27T 


Sk{Kx,Kv) 


■si 

SI 


{k  +  l/2)AKxtx 


(/  +  1/2)A  Kyty 
2 


—  erf 


—  erf 


(fc  -  1/2)  A  KJZ 


x 


(/  -  1/2)A Kyiv 
2 
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and  where  erf(-)  is  the  error  function.  The  accuracy  of  this  approximation  and  that  of 
a  few  similar  approximations  to  Equation  3.3  is  discussed  in  Appendix  B.  It  is  shown 
in  Appendix  B  that  Equation  3.5  is  generally  accurate  to  within  a  small  fraction  of 
a  percent. 


3.1.2  Random  Realizations  of  the  Impulse  Response  Function. 

The  impulse  response  functions  will  be  generated  one  delay  at  a  time  starting 
with  ihe  smallest  delay  (j  =  0)  and  working  to  the  largest  delay  (j  =  jmax)-  For  each 
delay,  the  steps  that  are  performed  are  outlined  in  this  subsection. 

The  angular  spectrum  of  the  impulse  response  function  in  the  jth  delay  bin 
at  the  output  of  the  mth  antenna  is 


hm{kAKz,lAKv,jAr) 


if  Eqn.  3.2 

0  otherwise 


(3.7) 


The  normalization  4n2 / AKxAKvAt  is  chosen  so  that  when  the  Fourier  transforms 
from  angle  to  space  are  complete,  the  quantity  hm(jAT,kAt)AT  will  represent  the  re¬ 
ceived  signal  (without  modulation  included)  during  the  delay  interval  jAr  to 
(j  +  l)Ar. 
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The  quantity  £(k,  l)  is  a  complex,  zero  mean,  Gaussian  random  variable 
with  the  properties 


«(M)>  =  0 

(£(MKK0))  =  o  (3.8) 


where  dkti  is  the  Kronecker  delta  symbol.  The  random  numbers  £  may  be  generated 
using  the  equation 


£  =  \J~lnUi  exp(27riu2)  (3.9) 

where  U!  and  u2  are  independent  random  numbers  uniformly  distributed  on  the  in- 
terval  [0,1). 

There  is  no  explicit  y  coordinate  since  realizations  will  be  generated  only  at  y 
coordinates  corresponding  to  antenna  centers.  This  is  accomplished  by  first  generating 
Nv  samples  of  hm  in  the  Kv  domain  and  then  by  Fourier  transforming  hm  to  the  y 
domain  using  discrete  Fourier  transforms  (DFTs).  One  DFT  is  performed  at  each 
Kx  coordinate  grid  point  for  each  unique  antenna  center  y  coordinate.  In  continuous 
notation,  this  transform  is 


~  r00  ~  dK 

hm(Kx,r)=  /  exp{iKvym)hm(Kx,Kv,T)— v-  (3.10) 

J— oo  Z7T 

where  ym  is  the  y  coordinate  of  the  mth  antenna  center.  The  DFT  equivalent  of 
Equation  3.10  is 


A  K  N^2 

hm(kAKx,jAr)  =  —^  £  exP(UAKvym)hm(kAKxtlAKvJAr)  .  (3.11) 

/7F  l=-N,/2 

Now  the  Kz  to  x  Fourier  transform  is  performed  for  each  antenna.  In  this 
case,  a  fast  Fourier  transform  (FFT)  may  be  used.  In  continuous  notation 
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The  discrete  equivalent  of  Equation  3.12  is 


A  K  N*!2~y 

hm(kAx,jAT )  =  — — -  53  exp(tV  AKxkAx)hm{k!  AKXJ  Ar)  (3.13) 

27r  k'=-Nx/2 

where  Nx  is  the  number  of  x  coordinate  grid  points. 

The  quantity  hm(kAx,jAr)  is  the  impulse  response  function  at  the  output  of 
the  mtk  antenna  when  that  antenna  is  located  at  x  position  kAx.  Under  the  frozen-in 
model,  a  time  step  At  may  be  associated  with  Ax  using  the  effective  velocity: 


At  =  Ax{t0H%)  .  (3.14) 

Thus  hm(kAx,jAT)  also  represents  the  impulse  response  function  at  the  output  of 
the  mth  antenna  at  time  kAt. 


3.1.3  Grid  Sizes. 

The  angular  grid  sizes  AKZ  and  AKV  are  chosen  so  the  magnitude  of  the 
antenna  filtered  GPSD  will  be  negligible  at  the  edges  of  the  angular  grid,  and  so  there 
will  be  a  sufficient  number  of  samples  per  decorrelation  distance  to  resolve  the  fades 
of  the  Rayleigh  amplitude  distribution. 

The  delay  grid  size  At  is  generally  chosen  on  the  basis  of  the  sampling 
period  of  the  receiver.  The  number  of  delay  samples  N&  of  the  realizations  is  then 
chosen  so  the  total  delay  spread  of  the  realization  NpAr  will  encompass  at  least  95 
percent  of  the  delayed  signal  energy. 

Ks  Grid 


The  x  grid  spacing  is 


Ax  —  t\x  /  No 


(3.15) 


where  t\x  is  the  minimum  (of  the  values  for  all  antennas)  filtered  x  direction  decorre¬ 
lation  distance.  The  quantity  No  is  the  number  of  samples  per  decorrelation  distance 
or  time  where  N0  should  have  a  value  of  10  or  larger  in  order  for  the  realization  to  ac¬ 
curately  represent  the  duration  and  separation  of  the  fades  of  the  Rayleigh  amplitude 
distribution  (Dana  1982,  1988). 

The  number  of  x  grid  points  Nx  is  subject  to  two  constraints:  First,  Nx 
should  be  a  power  of  2  so  a  fast  Fourier  transform  can  be  used  to  transform  from 
the  Kx  domain  to  the  x  domain.  Second,  Nx  should  be  greater  than  or  equal  to 
100  N0.  Equivalently,  there  should  be  at  least  100  x  direction  decorrelation  distances 
and  therefore  100  decorrelation  times  in  the  realization  in  order  for  the  realization  to 
represent  a  reasonable  sample  of  the  random  process. 

Here  reasonable  is  defined  in  terms  of  the  application.  Typically,  the  real¬ 
izations  are  used  to  exercise  simulations  of  trans-ionospheric  communications  links. 
Each  realization  with  Nx/N0  approximately  equal  to  100  can  be  used  to  calculate 
link  performance  averaged  over  100  decorrelation  times.  The  resulting  receiver  per¬ 
formance  measures  will  have  some  statistical  variation  due  to  the  finite  time  duration 
of  the  simulation.  This  variation  can  be  reduced  by  either  using  longer  realizations 
or  by  using  several  realizations  with  the  same  signal  parameters  but  with  different 
random  numbers  and  by  averaging  the  results  from  the  multiple  realizations.  This 
latter  approach  allows  an  estimat. ;  of  the  statistical  variation  in  the  performance  to 
be  made. 


Returning  to  the  grid  spacing,  the  x  direction  distance  spanned  by  the  re¬ 
alization  is  Lx  =  Nx  Ax  and  A Kx  is 


A  Kx  =  ~  = 


2n 


Nx  Ax 


(3.16) 


The  functional  dependence  of  the  antenna  filtered  GPSD  on  Kx  is  given  approximately 
by  exp[— {KxIax/2)2}.  In  discrete  notation,  Kx  =  kAKx  where  |&|  has  a  maximum 
value  of  Nx/2.  Using  the  maximum  value  for  Kx  of  NxAKx/2.,  the  value  of  the  GPSD 
at  the  edge  of  the  Kx  grid  is  proportional  to  exp[-(7riVo/2)2j  which  is  essentially  equal 
to  zero  when  Nq  equals  10. 
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Kv  Grid 

The  number  of  Kv  grid  points  is  determined  by  requiring  that  the  magni¬ 
tude  of  the  antenna  filtered  GPSD  is  small  at  the  edge  of  the  Kv  grid.  The  func¬ 
tional  dependence  of  the  antenna  filtered  GPSD  on  Kv  is  given  approximately  by 
exp(— (Ky^/2)2].  A  convenient  minimum  value  of  the  antenna  filtered  GPSD  at  the 
edge  of  the  Kv  grid  is  exp(-7r2).  Using  this  value  results  in  the  following  condition 
on  Nv: 


N  > 

*  -  AK/a, 


(3.17) 


where  iAv  is  the  minimum  (of  the  values  for  all  antennas)  filtered  y  direction  decor¬ 
relation  distance.  It  is  desirable  to  minimize  the  number  of  Kv  samples  in  order  to 
minimize  computer  memory  and  execution  time  requirements.  A  reasonable  minimum 
value  for  Ny  is  32. 

The  A Kv  grid  spacing  is  then 


A  K„  = 


27T 

Ty 


(3.18) 


where  the  length  Lv  is  chosen  to  be 


Ly  —  maximum  [16£/V,4maximum(|ym|)]  .  (3.19) 

The  first  condition,  Lv  =  16^,  results  from  letting  Nv  =  32  in  Equation  3.17.  The 
second  condition,  Lv  =  4  maximum  (|ym|),  is  chosen  to  minimize  aliasing  of  the  impulse 
response  functions  of  the  two  antennas  with  the  largest  separation  in  the  y  direction. 
If  the  value  of  Lv  is  determined  by  tins  latter  condition,  then  the  value  of  Ny  will  need 
to  be  larger  than  32  in  order  to  satisfy  the  condition  on  Nv  given  by  Equation  3.17. 

Delay  Grid 

The  delay  grid  size  is  usually  chosen  so  there  are  an  integer  number  of 
delay  samples  in  a  modulation  period.  The  number  of  delay  samples  is  chosen  so 
the  delay  spread  of  the  realization,  No  At,  will  encompass  at  least  95  percent  of  the 
delayed  energy.  Without  an  antenna,  95  percent  of  the  signal  energy,  under  isotropic 
scattering  conditions,  arrives  with  delays  in  the  range  -0.25  <  wC(J/,r  <  3.45  when  the 
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delay  parameter  a  is  set  at  its  minimum  value  of  4  (Wittwer  1980).  To  be  conservative, 
this  same  criterion  is  applied  under  all  conditions.  This  will  result  in  a  delay  spread 
that  encompasses  more  than  95  percent  of  the  energy  for  anisotropic  scattering  cases 
with  antennas.  The  minimum  number  of  delay  samples  required  is  then 


Nd  > 


3.7 

27r/x,mmAr 


(3.20) 


where  /a, min  is  the  minimum  (of  the  values  for  all  antennas)  filtered  value  of  the 
frequency  selective  bandwidth. 


3.2  TURBULENT  MODEL. 


The  starting  point  for  the  turbulent  model  is  the  angular  spectrum  of  the 
impulse  response  function.  At  discrete  time  nAt  (At  —  tq/No),  the  angular  compo¬ 
nents  of  the  impulse  response  function  at  the  output  of  the  mth  antenna  for  jth  delay 
bin  are 


hm(kAKz,lAKy,jAr,nAt) 


0  otherwise  . 

(3.21) 


The  difference  between  this  equation  and  Equation  3.7  is  that  the  random  numbers 
£(/:,/,  n)  are  now  functions  of  time  as  well  as  angle.  The  impulse  response  functions  at 
the  outputs  of  multiple  antennas  are  then  generated  from  Jim(kAKx,lAKv,jAT,nAt) 
by  performing  the  following  two  discrete  Fourier  transforms: 


A k  A  K 

hm(j Ar,nAt)  =  — — -  ]T  exp(ikAKxxm)  x 

47r  k=-Nx/2 

Ny/ 2-1 

exp(tlAKvym)hm(kAKx,lAKv,j  At  ,  nAt)  (3.22) 

l=-N,/2 
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where  xm  and  ym  are  the  x-y  coordinates  of  the  mth  antenna  position.  The  angular 
grid  size  A  Kz  for  this  model  is  chosen  using  the  same  criteria  that  were  used  above 
to  select  Ky. 

At  the  next  time  step,  (n  +  l)Af,  the  random  numbers  £(k,l,n)  are  recal¬ 
culated  according  to  the  Doppler  frequency  spectrum.  The  f~ 4  form  for  the  Doppler 
spectrum  is  convenient  because  the  random  numbers  £(k,l,n)  can  be  readily  synthe¬ 
sized  using  two-pole  low-pass  filters  (see,  for  example,  Wittwer  1980). 

Once  the  noise  samples  £  have  been  updated,  the  impulse  response  functions 
at  time  (n  +  1)  At  are  again  calculated  using  Equation  3.22. 
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SECTION  4 


EXAMPLES 


Several  examples  of  the  received  voltage  out  of  a  matched  filter  are  presented 
in  this  section.  These  examples  are  intended  to  illustrate  the  effects  of  frequency 
and  spatial  selectivity  and  the  differences  between  frozen-in  and  turbulent  models. 
Examples  for  specific  system  applications  may  be  found  in  Bogusch,  et.  al.  (1981) 
and  in  Bogusch,  Guigliano,  and  Knepp  (1983). 


4.1  MATCHED  FILTER  OUTPUT. 


The  output  of  a  matched  filter  can  be  constructed  by  convolving  the  im¬ 
pulse  response  function  of  the  channel  and  antenna  with  the  combined  impulse  re¬ 
sponse  function  of  the  transmitter  and  receiver.  A  second  approach  is  to  construct 
the  combined  frequency  response  of  the  transmitter,  channel,  antenna,  and  receiver 
and  then  to  Fourier  transform  that  result  to  obtain  the  matched  filter  output.  This 
latter  approach  is  used  for  the  examples  presented  in  this  section. 

The  starting  point  is  to  calculate  the  channel/antenna  transfer  function 
which  is  the  Fourier  transform  of  the  impulse  response  function: 

H(u,t)=  f  h(r,t)  exp(-iur)dT  (4.1) 

Jo 

This  function  represents  the  response  of  the  channel  and  antenna  at  time  f  to  a 
transmitted  sinewave  with  radian  frequency  u>. 

lor  a  transmitted  square  pulse  with  a  chip  duration  Tc,  the  voltage  out  of 
the  matched  filter  at  time  t  can  be  written  as 


/OO  (hj 

M{w)H(ut)  exp(iwr)—  (4.2) 

-oo  2tt 

where  r  is  the  time  delay  of  the  matched  filter  relative  to  the  nominal  time  of  arrival 
(e.g.  the  time  of  arrival  under  benign  propagation  conditions). 
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The  combined  spectrum  of  the  transmitted  square  pulse  and  the  receiver 
matched  filter  is 


M(w)  =  Tc 


sin2  (uTc  / 2) 

(wTc/ 2)J 


(4.3) 


The  impulse  response  function  is  generated  with  delay  samples  of  size 
A  r,  so  the  discrete  chr  .el/antenna  transfer  function  has  an  unambiguous  frequency 
response  over  2ir/Ar  radians.  If  this  bandwidth  is  divided  into  NF  frequency  samples, 
then  the  channel/ antenna  transfer  function  at  discrete  time  nAt  is 


nd- i 

II(kAu,  nAt)  =  [^(iAr,nAt)Ar)  exp  [— t(k  -  ATj? /2)Awj  Ar]  (4.4) 

j=o 

where  Aw  =  2-n/(NFAr).  The  range  of  the  index  k  in  this  equation  is  from  0  to 
NF  -  1  representing  a  range  of  frequencies  from  -NF Aw/2  to  ( NF  -  l)Aw/2.  Of 
course  the  number  of  frequency  samples  should  be  at  least  as  large  as  the  number 
of  delay  samples.  However,  if  the  number  of  delay  samples  was  chosen  to  be  the 
minimum  number  required  by  Equation  3.21,  then  it  may  be  necessary  to  chose  the 
number  of  frequency  samples  to  be  larger  that  the  number  of  delay  samples  in  order 
to  minimize  aliasing  (in  delay)  of  the  matched  filter  output. 

The  output  voltage  of  the  matched  filter  as  a  function  of  time  and  relative 
delay  is  given  by 


u(r,nAt) 


Aw  Tc  Np^~l  sin2  (k  AwTc  /  2) 
2?r  k=-NF/2  {kAu>Tc /2)2 


H (kAu>,  nAt)  exp(ikAuT) 


(4.5) 


If  At  is  chosen  to  be  Tc/ 2,  then  AwTc/2  =  2n/NF  and  u(r,t)  represents 
a  signal  that  is  band-limited  to  the  frequency  range  -l/Tc  to  +l/Tc.  Note  that 
the  matched  filter  output  u(r,t)  is  unambiguous  in  delay  over  the  interval  from  0  to 
(Nf  -  l)Ar  compared  to  the  delay  interval  of  0  to  (Np- 1)  for  the  original  ACIRF 
realization. 
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4.2  SCINTILLATION  EFFECTS  ON  MATCHED  FILTER  OUTPUT. 


In  the  examples  that  follow,  the  chip  rate  Rc  is  set  at  1  Mhz,  and  the  random 
realizations  of  the  impulse  response  function  are  generated  with  a  delay  sample  size  of 
Tc/ 2  (7c  =  1  /Rc)-  However,  the  frequency  selective  effects  depend  only  on  the  ratio 
of  the  frequency  selective  bandwidth  to  the  chip  rate  fo/Rc •  For  the  antenna  examples, 
a  uniformly  weighted  circular  antenna  and  isotropic  scattering  are  assumed.  Antenna 
effects  then  depend  only  on  the  ratio  of  the  antenna  diameter  D  to  the  decorrelation 
distance  Iq  and  the  antenna  pointing  angle. 


4.2.1  Frequency  Selective  Effects. 

In  a  high  data  rate  communications  link,  the  major  effect  of  frequency  se¬ 
lective  fading  is  intersymbol  interference.  Even  relatively  small  amounts  of  delay 
spread  can  catastrophically  degrade  demodulation  performance  in  such  a  link  using 
conventional  matched  filter  detection  techniques. 

Figure  4.1  shows  examples  of  the  matched  filter  output  amplitude  for  three 
levels  of  frequency  selective  propagation  disturbances,  characterized  by  the  ratio  of 
the  frequency  selective  bandwidth  /0  to  the  chip  rate  Rc. 

The  impulse  response  functions  were  generated  using  the  frozen-in  model 
and  a  small  antenna  (i.e.  D  <  to).  Each  frame  in  the  figure  provides  a  three  dimen¬ 
sional  picture  of  the  matched  filter  output  for  a  single  transmitted  pulse  as  a  function 
of  time  delay  (abscissa)  and  time  (scale  directed  into  the  figure).  The  total  duration 
of  each  of  the  frames  is  10  decorrelation  times. 

In  the  top  frame  the  frequency  selective  bandwidth  is  equal  to  the  chip  rate 
and  only  a  small  amount  of  distortion  is  evident  in  the  waveform  (which  is  slightly 
rounded  due  to  band  limiting  at  the  first  nulls  of  the  signal  spectrum).  The  effect  of 
fading  can  be  seen  in  this  frame  as  the  peak  amplitude  rises  and  falls  with  time.  Some 
minor  distortion  of  the  output  amplitude  is  seen  but  for  the  most  part  the  signal  is 
contained  within  the  period  of  one  chip.  This  channel  is  nearly  flat  fading  which  means 
that  all  frequency  components  within  the  signal  bandwidth  propagate  essentially  the 
same  way  through  the  disturbed  ionosphere.  There  is  very  little  time  delay  spread 
beyond  one  chip  in  the  matched  filter  output. 

The  middle  frame  in  Figure  4.1  shows  the  matched  filter  output  amplitude 
for  the  case  were  fo/Rc  is  equal  to  0.2.  For  this  smaller  value  of  the  frequency  selective 
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Figure  4.1.  Effects  of  frequency  selective  fading. 
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bandwidth,  more  of  the  signal  energy  is  arriving  with  delays  of  more  than  one  chip, 
and  there  are  multiple  distinct  peaks  in  the  matched  filter  output  amplitude.  It  is 
these  structures  that  can  cause  delay  tracking  algorithms  to  loose  lock  and  that  cause 
intersymbol  interference  which  can  degrade  demodulation  performance. 

The  bottom  frame  shows  a  highly  disturbed  case  where  f0  is  a  tenth  of 
the  chip  rate.  This  causes  signal  energy  to  be  spread  over  approximately  eight  chip 
periods.  When  a  contiguous  set  of  pulses  is  transmitted,  the  delay  spread  of  the 
received  signal  results  in  the  simultaneous  reception  of  information  from  about  eight 
previous  chips  which  can  produce  severe  intersymbol  interference. 

An  effect  due  to  the  frozen-in  model  that  is  evident  in  Figure  4.1  is  that 
the  signal  arriving  at  long  delays  varies  more  rapidly  in  time  than  the  signal  arriving 
at  shorter  delays.  A  comparison  of  the  matched  filter  output  amplitude  generated 
with  the  frozen-in  and  turbulent  models  is  shown  in  Figure  4.2  for  the  case  where 
fo/Rc  is  equal  to  0.1.  The  top  frame  in  this  figure  for  the  frozen-in  model  is  just 
a  reproduction  of  the  bottom  frame  in  Figure  4.1.  Again  10  decorrelation  times  of 
the  signal  are  plotted.  It  should  be  recalled  when  comparing  the  two  frames  that 
the  Doppler  spectrum  for  the  frozen-in  model  has  a  Gaussian  form,  and  the  Doppler 
spectrum  for  the  turbulent  model  has  the  /“4  form.  It  is  this  fact  that  accounts  for 
the  more  spiky  appearance  of  the  turbulent  model  output  amplitude.  The  important 
difference  between  the  two  frames  is  that  the  turbulent  model  amplitude  has  the  same 
fading  rate  at  all  delays. 


4.2.2  Spatially  Selective  Effects. 

Spatially  selective  effects  are  important  for  high  data  rate  communications 
links  that  rely  on  large  antennas  to  achieve  sufficient  signal-to-noise  ratios  for  low 
error  rate  data  demodulation.  Scattering  loss  of  an  antenna  is  a  function  of  the  size  of 
the  antenna  D  relative  to  the  decorrelation  distance,  to.  When  to  is  greater  than  D, 
the  electric  field  is  highly  correlated  across  the  face  of  the  antenna  and  the  full  gain 
of  the  antenna  is  realized.  However,  the  antenna  may  be  located  at  a  position  where 
the  incident  power  is  in  a  deep  fade.  The  solution  to  this  problem  is  to  have  multiple 
antennas  physically  separated  by  a  distance  larger  than  the  maximum  decorrelation 
distance.  The  probability  of  having  all  antennas  simultaneously  experience  deep  fades 
in  the  received  power  is  then  substantially  reduced. 

The  problem  of  spatial  selectivity  occurs  when  to  is  less  than  D  and  the 
electric  field  is  decorrelated  across  the  face  of  the  antenna.  In  this  case,  the  induced 
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Figure  4.2.  Comparison  of  matched  filter  output  amplitude  for  the 
frozen-in  and  turbulent  models. 
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voltages  in  the  antenna  add  noncoherently  due  to  the  random  phase  variations  in 
the  electric  field,  and  a  loss  of  signal  power,  or  equivalently  of  antenna  gain,  is  the 
result.  From  another  perspective,  this  loss  occurs  when  the  angular  scattering  process 
responsible  for  amplitude  and  phase  scintillation  and  frequency  selective  effects  also 
causes  some  of  the  transmitted  signal  energy  to  be  scattered  out  of  the  antenna  beam. 

Figure  4.3  shows  examples  of  the  matched  filter  output  amplitude  for  three 
levels  of  spatially  selective  propagation  disturbance,  characterized  by  the  ratio  of  the 
antenna  size  to  the  decorrelation  distance.  The  ratio  of  the  frequency  selective  band¬ 
width  to  the  chip  rate  is  0.1,  and  the  antenna  is  pointing  along  the  line  of  sight.  The 
top  frame  is  for  the  case  where  Iq  is  much  greater  than  D,  and  is  just  a  reproduction  of 
the  bottom  frame  of  Figure  4.1.  The  middle  frame  is  for  a  lo/D  ratio  of  0.5  where  the 
scattering  loss  is  3.1  dB.  The  effect  of  the  antenna  is  to  preferentially  attenuate  the 
signal  energy  arriving  at  large  angles  and  also  at  large  delays  and  thereby  to  reduce 
the  delay  spread  of  the  output  signal.  In  the  bottom  frame  where  lo/D  is  equal  to  0.2, 
the  output  signal  is  almost  flat  with  very  little  delay  spread  distortion  of  the  matched 
filter  output.  Although  this  substantially  reduces  the  effects  of  frequency  selective 
fading,  the  cost  is  an  8.8  dB  reduction  in  the  average  signal  power. 

Finally,  Figure  4.4  shows  examples  of  the  matched  filter  output  amplitude 
for  three  values  of  the  pointing  angle  ip.  The  ratio  of  the  frequency  selective  bandwidth 
to  the  chip  rate  is  0.1,  and  the  ratio  of  decorrelation  distance  to  antenna  diameter  is 
0.5.  The  top  frame  for  a  pointing  angle  of  zero  is  just  a  reproduction  of  the  middle 
frame  of  Figure  4.3.  The  average  scattering  loss  for  this  case  is  3.1  dB.  The  bottom 
two  frames  show  the  matched  filter  output  amplitude  for  pointing  angles  of  one-half 
beamwidth  (ip  =  60/ 2)  with  a  scattering  loss  of  4.6  dB  and  one  beamwidth  (ip  =  0O) 
with  a  scattering  loss  of  9.2  dB. 

Although  the  average  scattering  losses  are  about  the  same,  the  bottom  frame 
of  Figure  4.4  (lo/D  =  0.5  and  ip  =  60 )  and  the  bottom  frame  of  Figure  4.3  (lo/D  =  0.2 
and  ip  =  0)  are  qualitatively  quite  different.  For  the  case  with  the  pointing  angle  equal 
to  a  beamwidth,  the  received  power  is  much  more  spread  out  in  delay  compared  to 
the  case  with  zero  pointing  angle  where  the  signal  energy  is  concentrated  near  zero 
delay.  This  is  due  to  the  fact  that  the  antenna  pointed  away  from  the  line  of  sight  has 
relatively  higher  gain  at  large  angles  and  long  delays  and  relatively  lower  gain  at  small 
angles  and  short  delays  than  does  an  antenna  pointed  along  the  line  of  sight.  Thus 
for  an  antenna  pointed  away  from  the  line  of  sight,  increased  scattering  loss  does  not 
necessarily  result  in  reduced  frequency  selective  effects. 
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Figure  4.3.  Effects  of  spatially  selective  fading. 
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Figure  4.4.  Effects  of  beam  pointing. 
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APPENDIX  A 


INPUTS  AND  OUTPUTS  OF  ACIRF 


A.l  ACIRF  INPUT  FILE. 


The  key  inputs  to  ACIRF  are  listed  in  Table  A.l.  These  include  channel 
parameters  which  describe  the  second  order  statistics  of  the  signal  incident  on  the  face 
of  any  antenna,  antenna  parameters,  and  realization  parameters. 

There  are  five  channel  models  included  in  ACIRF:  frozen-in  model  with 
antennas  (IFADE  =  1);  turbulent  model  with  antennas  (IFADE  =  2);  frozen-in  model 
without  antennas  (IFADE  =  3);  turbulent  model  without  antennas  (IFADE  -  4);  and 
a  flat  fading  model  (IFADE  =  5).  Models  3  and  5  are  also  implemented  in  Wittwer’s 
channel  model  CIRF  (Wittwer  1980). 

The  key  realization  parameter  inputs  are  the  number  of  delay  samples  and 
the  delay  sample  size.  The  delay  sample  size  should  be  no  larger  than  one-half  of 
the  chip  period.  For  most  applications,  setting  A r  equal  to  Tc/ 2  is  sufficient.  The 
number  of  delay  samples  required  is  determined  from  Equation  3.21. 

An  example  input  file  is  listed  in  Table  A.2.  This  input  file  was  used  to 
generate  realizations  of  the  impulse  response  function  from  which  the  plots  in  Figure 
4.4  were  created.  Note  that  comment  lines  in  the  input  file  start  with  a  semicolon  (;). 


A. 2  FORMATTED  ACIRF  OUTPUT  FILE. 


The  formatted  output  file  generated  using  the  example  ACIRF  data  file  is 
listed  in  Table  A. 3a.  The  first  page  of  this  file  provides  a  summary  of  the  inputs  and 
lists,  for  each  antenna,  ensemble  values  of  scattering  loss,  and  filtered  frequency  se¬ 
lective  bandwidth,  decorrelation  distances,  and  decorrelation  time.  Subsequent  pages 
of  the  output  file  list,  for  each  antenna,  the  measured  values  of  these  quantities.  The 
last  page  lists  ensemble  and  measured  values  for  the  cross  correlation  of  the  signal  out 
of  the  antennas.  The  algorithms  used  to  measure  these  quantities  are  given  in  Dana 
(1986). 
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Table  A.l.  Key  ACIRF  inputs. 


Channel  Parameters: 

Channel  model,  IFADE 

Frequency  Selective  Bandwidth,  /o 

(Hertz) 

Decorrelation  Distances,  iz  and  ly 

(meters) 

Decorrelation  Time,  tq 

(seconds) 

Antenna  Parameters: 

Beamwidths,  0OU  and  6ov 

(degrees) 

Positions,  u  and  v 

(meters) 

Rotation  angle,  x 

(degrees) 

Pointing  angle,  ip 

(degrees) 

Pointing  azimuth,  £ 

(degrees) 

Realization  Parameters: 

Number  of  delay  samples,  ND 

Delay  sample  size,  At 

(seconds) 

For  each  delay  bin,  amplitude  moments  and  the  scintillation  index  S4  of 
the  impulse  response  function  are  measured  and  output  in  the  formatted  file.  The 
ensemble  power  in  the  jth  delay  bin  for  the  mth  antenna  is 


Pmj  ~ 


/•(;+l)Ar 

j'Ar 


Pm[r)dT 


for  j  =  0,  1,  •  •  •,  Nd~  1.  The  power  impulse  response  function  Pm(r)  is 


(A.1) 


•«=/:£/: 


r°°  dKu 


27T 


OJK.  -  K0„K,  -  K.,)Sk{K',K„T) 


(A.2) 


and  is  given  in  Frasier  (1989).  The  ensemble  values  of  the  amplitude  moments  are 
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Table  A.2.  Example  ACIRF  input  file 


THIS  DATA  FILE  IS  FOR  ACIRF  USER'S  GUIDE  EXAMPLES 

******  CASE  NUMBER 

KASE 

1001/ 

******  ALPHANUMERIC  IDENTIFICATION  (80  CHARACTERS  OR  LESS  ENCLOSED  IN  '  ') 

ACIRF  USERS  GUIDE  EXAMPLES' 

******  CHANNEL  PARAMETERS 

IFADE  FO(HZ)  TAUO(S)  DLX(M)  DLY(M) 

1  0.1E6  3.0E-3  5.0  5.0/ 

******  ANTENNA  PARAMETERS  (READ  ONLY  IF  IFADE  -  1  OR  IFADE  *  2) 

NUMANT  FREQC(HZ) 

3  2.99793E9/ 

BEAMWIDTHS.  POSITIONS,  ROTATION  ANGLES,  AND  POINTING  ANGLES 


(ONE  LINE  FOR  EACH  ANTENNA) 

;  BWU(DEG) 

BWV(DEG) 

UPOS(M) 

VPOS(M) 

CHI (DEG) 

PSI(DEG) 

EPS (DEG) 

0.5890 

0.5896 

0.0 

0.0 

00.0 

0.0 

0.0/ 

0.6896 

0.5896 

0.0 

0.0 

00.0 

0.2948 

0.0/ 

0.5896 

0.5896 

0.0 

0.0 

00.0 

0.5896 

0.0/ 

;******  REALIZATION  PARAMETERS 

NDELAY  DELTAU(S) 

NTIMES 

NKX 

NKY 

ISEED 

NO 

16 

5.0E-7 

1024 

1024 

33 

9771975 

10/ 

-3 


Table  A. 3a.  Example  ACIRF  output  file  (page  1,  summary  of  input  and 
ensemble  values). 

ACIRF  CHANNEL  SIMULATION 
CASE  NUMBER  1001 

TEMPORAL  VARIATION  FROM  FROZEN- IN  MODEL 
REALIZATION  IDENTIFICATION: 

ACIRF  USERS  GUIDE  EXAMPLES 

CHANNEL  PARAMETERS 

FREQoENCT  SELECTIVE  BANDWIDTH  (HZ)  ■  1.000E+06 

DECORRELATION  TIME  (SEC)  -  3.000E-03 

X  DECORRELATION  DISTANCE  (M)  -  6. OOOE+OO 

Y  DECORRELATION  DISTANCE  (M)  -  6. OOOE+OO 


REALIZATION  PARAMETERS 

NUMBER  OF  DELAY  SAMPLES  -  18 
DELAY  SAMPLE  SIZE  (SEC)  -  5.000E-07 
NUMBER  OF  TEMPORAL  SAMPLES  -  1024 
NUMBER  OF  KX  SAMPLES  -  1024 
NUMBER  OF  KY  SAMPLES  •  33 
INITIAL  RANDOM  NUMBER  SEED  -  0771075 
NUMBER  OF  TEMPORAL  SAMPLES  PER  TAUO  -  10 


ANTENNA  PARAMETERS 

NUMBER  OF  ANTENNAS  -  3 

CARRIER  FREQUENCY  (HZ)  -  2.008E+00 

ANTENNA  BEAMWIDTHS,  ROTATION  ANGLES  AND  POINTING  ANGLES 

N  BWU(DEG)  BWV(DEG)  CHI (DEG)  PSI(DEG)  EPS (DEG) 

1  .500  .500  .000  .000  .000 

2  .500  .500  .000  .206  .000 

3  .500  .500  .000  .500  .000 

ANTENNA  POSITIONS  IN  U-V  AND  X-Y  COORDINATES 

N  UPOS(M)  VPOS(M)  XPOS(M)  YPOS(M) 

1  OOOE+OO  .OOOE+OO  .OOOE+OO  .OOOE+OO 

2  .OOOE+OO  .OOOE+OO  .OOOE+OO  .OOOE+OO 

3  .OOOE+OO  .OOOE+OO  .OOOE+OO  .OOOE+OO 

ENSEMBLE  CHANNEL  PARAMETERS  AT  ANTENNA  OUTPUTS 

N  LOSS(DB)  POWER  FA(HZ)  ALX(M)  ALY(M)  TAUA(SEC) 

1  3.141  .485187  2.081E+06  7.178E+00  7.178E+00  4.307E-03 

2  4.802  .346811  1.574E+06  7.178E+00  7.178E+00  4.307E-03 

3  8.983  .128390  1.060E+06  7.178E+00  7.178E+00  4.307E-03 


Table  A.3b.  Example  ACIRF  output  file  (page  2,  summary  for 
antenna  1). 


MEASURED  PARAMETERS  FOR  REALIZATION/ ANTENNA  1 

MOMENTS  OF  VOLTAGE  AMPLITUDE  VERSUS  DELAY 
NORMALIZED  TO  ENSEMBLE  VALUES 
POW(J)  *  ENSEMBLE  POWER  IN  J-TH  DELAY  BIN 
J  »  T:  STATISTICS  OF  COMPOSITE  SIGNAL 


J 

POW(J) 

<A> 

<A**2> 

<A**3> 

<A**4> 

S4 

0 

2.313E-01 

.0668 

.0600 

.0780 

1.0344 

1 . 1368 

1 

1.210E-01 

1.0182 

.0707 

.8064 

.7833 

.7061 

2 

6.334E-02 

1.0642 

1.1008 

1 . 1242 

1.1372 

.0366 

3 

3.316E-02 

.0883 

.0866 

.0018 

.0032 

1.0203 

4 

1 . 736E-02 

.0261 

.8661 

.8137 

.7676 

1.0262 

S 

0.078E-03 

.0740 

.0333 

.8820 

.8206 

.0404 

6 

4.761E-03 

.0002 

.0610 

.8017 

.8188 

.0006 

7 

2.486E-03 

.0874 

1.0000 

1.0246 

1.0626 

1.0404 

8 

1.301E-03 

.0603 

.0303 

.0276 

.0160 

1.0373 

0 

6.810E-04 

.8717 

.7371 

.6008 

.4062 

.0071 

10 

3 . 564E-04 

1.0041 

1.0300 

1.0711 

1.1248 

1.0686 

11 

1 .866E-04 

1.0206 

1.0373 

1.0646 

1.0760 

.0001 

12 

0.761E-06 

.0761 

.0463 

.0068 

.8670 

.0660 

13 

K . 108E-06 

1 . 1066 

1.1000 

1.2640 

1.3378 

.0431 

14 

2 . 673E-06 

.0104 

.8846 

.8006 

.0607 

1.2063 

16 

1.300E-06 

1.0141 

1.0324 

1.0470 

1.0663 

.0010 

T 

4.352E-01 

.0073 

.0661 

.0071 

.8270 

.8808 

EASURED  REALIZATION  SIGNAL  PARAMETERS 

MEAN  POWER  OF  REALIZATION  -  .473600 
TOTAL  SCATTERING  LOSS  (DB)  •  3.24S 
FREQUENCY  SELECTIVE  BANDWIDTH  (HZ)  -  2.132E+06 
X  DIRECTION  DECORR.  DISTANCE  (M)  -  7.326E+00 
NUMBER  OF  SAMPLES  PER  DECORR.  DIST  -  10.206 
DECORRELATION  TIME  (SEC)  -  4.306E-03 

EAN  POWER  IN  GRID 

POWER  IN  KX-KY  GIRD  •  .484762 
POWER  LOSS  OF  GRID  (DB)  -  .004 
POWER  IN  DELAY  GRID  -  .486162 
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Table  A.3c.  Example  ACIRF  output  file  (page  3,  summary  for 
antenna  2). 


MEASURED  PARAMETERS  FOR  REALIZATION/ANTENNA  2 

MOMENTS  OF  VOLTAGE  AMPLITUDE  VERSUS  DELAY 
NORMALIZED  TO  ENSEMBLE  VALUES 
POW(J)  -  ENSEMBLE  POWER  IN  J-TH  DELAY  BIN 
J  -  T:  STATISTICS  OF  COMPOSITE  SIGNAL 


J 

POW(J) 

<A> 

<A**2> 

<A**3> 

<A**4> 

S4 

0 

1 . 280E-01 

‘.9974 

1.0276 

1.0774 

1 . 1387 

1.0766 

1 

8.267E-02 

1.0332 

1.0292 

.9922 

.0301 

.8696 

2 

6.231E-02 

1.0047 

1.0114 

1.0137 

1.0034 

.9807 

3 

3.267E-02 

.9067 

.8624 

.8634 

.8699 

1 . 1672 

4 

2.017E-02 

.9006 

.8167 

.7400 

.6690 

1.0029 

6 

1.233E-02 

.8829 

.7841 

.6963 

.6170 

1.0036 

6 

7.476E-03 

.9827 

.9663 

.9247 

.8834 

.9663 

7 

4.601E-03 

.0811 

.9739 

.9699 

.9628 

1.0160 

8 

2.603E-03 

1.0063 

1.0322 

1.0719 

1.1137 

1.0442 

9 

1.603E-03 

.7867 

.6194 

.4786 

.3686 

.9324 

10 

9.496E-04 

1.0011 

1.0166 

1.0284 

1.0236 

.9906 

11 

6.600E-04 

1.0708 

1.0918 

1.0803 

1.0432 

.8661 

12 

3 . 290E-04 

1.0146 

.9620 

.8434 

.7138 

.7683 

13 

1.926E-04 

1.2411 

1.4731 

1.6962 

1.9090 

.8716 

14 

1 . 124E-04 

1.0238 

.9946 

.9298 

.8416 

.8377 

16 

6.538E-06 

.9993 

.9698 

.8861 

.7892 

.8447 

T 

3 . 466E-01 

.9600 

.9090 

.8600 

.7840 

.9474 

IASURED  REALIZATION  SIGNAL  PARAMETERS 


MEAN  POWER  OF  REALIZATION  -  .341379 
TOTAL  SCATTERING  LOSS  (DB)  -  4.068 
FREQUENCY  SELECTIVE  BANDWIDTH  (HZ)  -  1.637E+06 
X  DIRECTION  DECORR.  DISTANCE  (M)  -  8.476E+00 
NUMBER  OF  SAMPLES  PER  DECORR.  DIST  -  11.808 
DECORRELATION  TIME  (SEC)  -  5.086E-03 

LAN  POWER  IN  GRID 


POWER  IN  KX-KT  GIRD 
POWER  LOSS  OF  GRID  (DB) 
POWER  IN  DELAY  GRID 


.346241 

.006 

.346621 


Table  A.3d.  Example  ACIRF  output  file  (page  4,  summary  for 
antenna  3). 


MEASURED  PARAMETERS  FOR  REALIZATION/ ANTENNA  3 

MOMENTS  OF  VOLTAGE  AMPLITUDE  VERSUS  DELAY 
NORMALIZED  TO  ENSEMBLE  VALUES 
POW(J)  -  ENSEMBLE  POWER  IN  J-TH  DELAY  BIN 
J  -  T:  STATISTICS  OF  COMPOSITE  SIGNAL 


J 

POW(J) 

<A> 

<A**2> 

<A**3> 

<A**4> 

S4 

0 

2.136E-02 

1.0277 

1.0802 

1 . 1302 

1.1686  1.0015 

1 

2. 181E-02 

1.0508 

1.0626 

1.0300 

.0706 

.8482 

2 

1.0B8E-02 

.0002 

.0000 

.0007 

.0674 

.0600 

3 

1.03OE-O2 

.8718 

.8152 

.8023 

.8213  1.2132 

4 

1.201E-02 

.8010 

.7008 

.0707 

.6000  1.0373 

5 

0.861E-03 

.8753 

.7767 

.6851 

.6062 

1.0058 

6 

7.308E-03 

.0670 

.0460 

.0260 

.8064  1.0038 

7 

5 . 300E-03 

.0001 

1.0160 

1.0123 

.0838 

.0530 

8 

3.772E-03 

1.0020 

1.0523 

1 . 1282 

1.2166  1.0034 

0 

2.C43E-03 

.7564 

.5634 

.4070 

.2825 

.8834 

10 

1.827E-03 

1.0246 

1.0423 

1.0436 

1.0173 

.0342 

11 

1 . 248E-03 

1.0000 

1 . 1060 

1.0637 

.0770 

.7730 

12 

8.430E-04 

1.0133 

.0423 

.8212 

.6776 

.7264 

13 

6.053E-O4 

1.2884 

1.5853 

1.8641 

2.1006 

.8230 

14 

3.760E-O4 

1.0026 

1.0404 

.0860 

.8020 

.7876 

16 

2.477E-04 

1.0060 

.0720 

.8080 

.7022 

.8228 

T 

1.260E-01 

.0246 

.8404 

.7772 

.7080 

.0820 

MEASURED  REALIZATION  SIGNAL  PARAMETERS 


MEAN  POWER  OF  REALIZATION 
TOTAL  SCATTERING  LOSS  (DB) 

FREQUENCY  SELECTIVE  BANDWIDTH  (HZ) 

I  DIRECTION  DECORR.  DISTANCE  (M) 
NUMBER  OF  SAMPLES  PER  DECORR.  DIST 
DECORRELATION  TIME  (SEC) 

MEAN  POWER  IN  GRID 

POWER  IN  KX-KY  GIRD  •  .125820 

POWER  LOSS  OF  GRID  (DB)  -  .010 

POWER  IN  DELAY  GRID  ■  .125033 


-  .110003 

-  0.212 

-  1.061E+06 

■  7 . 308E+00 

-  10.181 

■  4 . 386E-03 


Table  A.3e.  Example  ACIRF  output  file  (page  7,  ensemble  and  measured 
antenna  output  cross  correlation  coefficients). 


ENSEMBLE  ANTENNA  OUTPUT  CROSS  CORRELATION 

AMPLITUDE  OF  CROSS  CORRELATION 

N 

AMP(N-l) 

AMP(N-2) 

AMP(N*3) 

1 

1.000000 

.914081 

.699883 

2 

.914061 

1.000000 

.914666 

3 

. 099883 

.914065 

1.000000 

PHASE 

(RADIANS) 

OF  CROSS 

CORRELATION 

N 

PHS(N-l) 

PHS(N-2) 

PHS(N-3) 

1 

.000000 

.000000 

.000000 

2 

.000000 

.000000 

.000000 

3 

.000000 

.000000 

.000000 

MEAUSRED  ANTENNA 

OUTPUT  CROSS  CORRELATION 

AMPLITUDE  OF  CROSS  CORRELATION 

N 

AMP(N-l) 

AMP(N-2) 

AMP(N*3) 

1 

1.000000 

.910112 

. 703803 

2 

.910112 

1.000000 

.914799 

3 

. 703803 

.914799 

1.000000 

PHASE 

(RADIANS) 

OF  CROSS 

CORRELMION 

N 

PHS(N-l) 

PHS(N-2) 

PHS(N-3) 

1 

.000000 

.000000 

.000000 

2 

.000000 

.000000 

.000000 

3 

.000000 

.000000 

.000000 

FINAL  RANDOM  NUMBER  SEEDS 

ISEED 

m 

26650 

JSEED 

m 

63489 
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<  dj  > 

—  2  sJ^Pmj 

A 

V 

of 

II 

(A.3) 

<*?> 

=  j\/,rP-; 

<  a)  > 

=  2  Pi;  . 

mj 

The  mean  power  in  the  Kx  —  Ky  grid,  for  the  mth  antenna,  is 

Ni/2-l 

",/ 2-1 

PKG.m  =  E 

E  , 

(A.4) 

k=-N,/ 2 

l=~N,/  2 

and  the  mean  power  in  the  delay  gird,  for  the  mtk  antenna,  is 

ND- 1 

PDG.m  =  E  P*  •  (A.5) 

;=o 

These  last  two  quantities  are  computed  and  output  to  ensure  that  the  angular-delay 
grids  contain  most  of  the  signal  energy. 


A. 3  UNFORMATTED  ACIRF  OUTPUT  FILES. 


The  overall  file  structure  of  the  unformatted  ACIRF  output  files  which  con¬ 
tain  the  impulse  response  function  realizations  is  described  in  Table  A.4. 

The  first  record  of  an  unformatted  file  is  the  80  character  identification  that 
is  input.  The  Fortran  code  that  writes  this  record  is  shown  in  Table  A.5. 

The  problem  definition  data  record  A  contains  a  summary  of  the  input  data, 
and  Record  B  contains  a  detailed  definition  of  the  realizations.  The  generating  Fortran 
code  is  listed  in  Table  A.6.  A  description  of  the  important  words  in  Record  A  is  given  in 
Table  A. 7.  Because  a  common  output  file  format  is  used  by  several  channel  simulation 
programs  developed  at  Mission  Research  Corporation,  ACIRF  does  not  utilize  all  of 
the  words  in  the  record.  Only  the  words  that  contain  information  relevant  to  the  use 
of  ACIRF  realizations  are  listed  in  the  table. 
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Table  A.4.  Structure  of  unformatted  ACIRF  output  files. 


Description 

1 

Character  identification  record 

2 

Floating  point  record  A  (problem  definition  data) 

3 

Floating  point  record  B  (detailed  problem  definition  data) 

4 

Record  A 

5 

Floating  point  realization  data  record  Ci 

6 

Record  A  (if  necessary) 

7 

• 

Floating  point  realization  data  record  C2  (if  necessary) 

• 

• 

2n 

• 

Record  A  (if  necessary) 

2n+2 

Floating  Point  realization  data  record  C„  (if  necessary) 

Table  A. 5.  Fortran  code  that  generates  character  identification  record. 


PARAMETER  (NDENT=80) 

CHARACTER  IDENT*80 

WRITE (IUNIT)  NDENT , IDENT ( 1 : NDENT) 


Table  A. 6.  Fortran  code  that  generates  floating  point  header  records  A 
and  B. 


PARAMETER  (NDATA1=25 , NDATA2=32) 
C0MM0N/HEADER/RDATA1 (NDATAl) ,RDATA2(NDATA2) 
WRITE(IUNIT)  NDATAl , (RDATA1 (II) .11=1 .NDATAl) 
WRITE (IUNIT)  NDATA2 , (RDATA2 (II) ,11=1 .NDATA2) 
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Table  A. 7.  Description  of  header  Record  A. 


RDATA1(  1)  =  2.0 
RDATA1(  2)  =  KASE 
RDATA1(  3)  =  FREQC 
RDATA1(  4)  -  TAUO 
RDATA1(  5)  =  F0 
RDATA1(  6)  =  DLX 
RDATA1(  7)  =  DLY 
RDATA1(  9)  =  1.0 
RDATA1(13)  =  TDUR 
RDATA1(14)  =  NTIMES 
RDATAl(15)  =  DELTAT 
RDATA1(16)  =  NO 
RDATA1(20)  =  NDELAY 
RDATA1(22)  =  DELTAU 
RDATA1(23)  =  ISEED 
RDATAl(25)  =  MAXBUF 


ACIRF  Realization 
Case  Number 
Carrier  Frequency  (Hz) 

Decorrelation  Time,  tq  (sec) 

Frequency  Selecti-.e  Bandwidth,  fo  (Hz) 
x  Decorrelation  Distance,  tx  (m) 
y  Decorrelation  Distance,  lv  (m) 
Scintillation  Index,  S4 
Time  Duration  of  Realization,  (iVy  At) 
Number  of  Time  Samples,  Nt 
Time  Sample  Size,  At 
Samples  per  Decorrelation  Time,  N0 
Number  of  Delay  Samples,  No 
Delay  Sample  Size,  A r 
Random  Number  Seed 
Maximum  Buffer  Size 


The  ACIRF  subroutines  WRITER  and  READER  which  write  and  read 
the  unformatted  files  are  listed  in  Tables  A. 8  and  A.9,  respectively.  The  subroutine 
READER  may  be  adapted  to  read  the  impulse  response  function  into  a  link  simula¬ 
tion.  As  currently  written,  READER  outputs,  through  the  call  statement,  an  array 
containing  the  impulse  response  function  h(jAr,nAt)  ( j  =  1,  •  •  •  ,No)  at  a  fixed  time 
sample  nAt  (n  =  1,  •  •  • ,  iVr).  Consecutive  calls  of  READER  produce  consecutive  time 
samples  of  the  impulse  response  function. 

The  subroutine  WRITER  is  called  by  ACIRF  once  each  time  step  nAt  (n  = 
1  and  passed  an  array  of  the  impulse  response  function  h(jAr,nAt,)  ( j  = 

l,-1-,^).  WRITER  buffers  these  arrays  until  the  end  of  the  realization  is  reached 
or  until  the  number  of  impulse  response  function  arrays  equals  the  capacity  of  the 
buffer.  The  maximum  number  of  arrays  in  buffer  is  equal  to  the  largest  integer  which 
is  less  than  or  equal  to  (MAXBUF/2)/NDELAY.  The  ACIRF  parameter  MAXBUF  is 
the  maximum  number  of  real  words  in  an  unformatted  file  record.  The  current  value 
of  MAXBUF  is  4096.  Once  the  buffer  is  filled  to  capacity  or  the  end  of  the  realization 
is  reached,  WRITER  writes  the  buffer  into  the  file  as  a  single  record. 
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Table  A.8.  Listing  of  ACIRF  subroutine  WRITER 


SUBROUTINE  WRITER(RESET,NDEUT,NTIMES,NFILE,  IUNIT.H) 

C 

C* ******** ************  VERSION  1.0  ***  26  NOV  1987  *************************** 

C 

C  THIS  SUBROUTINE  WRITES  A  COMPLEX  ARRAT  INTO  A  FILE  ON  UNIT-IUNIT 

C 

C 

C  INPUTS: 

C 

C  RESET  -  LOGICAL  FUG  -  INTIALIZATION  IF  RESET  -  .TRUE. 

C  H  -  C0MPLEX*8  ARRAT 

C  HD  RUT  -  NUMBER  OF  DEUT  SAMPLES  IN  IMPULSE  RESPONSE  FUNCTION 

C  NTIMES  -  NUMBER  OF  TIME  SAMPLES  IN  IMPULSE  RESPONSE  FUNCTION 

C  NFILE  -  FILE  NUMBER 

C  I UNIT  ■  UNIT  NUMBER  OF  FILE 
C 

C  COMMON  BLOCK  INPUTS:. 

C 

C  RDATA1  -  HEADER  RECORD 
C 

C  PARAMETERS: 

C 

C  MAXBUF  -  MAXIMUM  BUFFER  SIZE  (RECORD  SIZE  IN  FILE  -  MAXBUF  +  1) 

C  MTENNA  •  MAXIMUM  NUMBER  OF  ANTENNAS  (EQUAL  TO  MAXIMUM  NUMBER  OF  FILES) 

C  NDATA1  -  SIZE  OF  HEADER  RECORD  1 

C  NDATA2  -  SIZE  OF  HEADER  RECORD  2 

C 

C***************************************************************************** 

c 

PARAMETER  (MAXBUF-409fl,MTENNA-8,NDATAl-25,NDATA2«32) 

LOGICAL  RESET 

COMMON/HEAD ER/RDATA1 (NDATA1) , RDATA2 (MTENNA , NDATA2) 

COMPLEX  H(NDEUT)  ,BUFFER(MAXBUF/2  .MTENNA) 

DIMENSION  NUMBER (MTENNA) 

C 

C******  IF  RESET  .EQ.  TRUE  THEN  INITIALIZE 
C 
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Table  A. 8.  Listing  of  ACIRF  subroutine  WRITER  (Concluded). 


IF (RESET) THEN 

NTBUF  -  (MAXBUF/2)/NDELAT 
NCXWRD  ■  NTBUF *NDELAT 
NRLWRD  -  2*NCXWRD 
DO  10  NN-l.MTENNA 
NUKBER(NN)  -  0 
10  CONTINUE 

RESET  -  .FALSE. 

RETURN 
END  IF 
C 

C  COUNT  NUMBER  OF  CALLS  -  ONE  PER  TIME  SAMPLE 
C 

}IW3ER(NFILE)  -  NUMBER(NFILE)  ♦  1 

JOFFST  -  MOD ( (NUMBER(NFILE) - 1 ) *NDELAT .NCXWRD) 

C 

C  WRITE  DELAY  SAMPLES  INTO  BUFFER 
C 

DO  20  JJ-l.NDELAY 

BUFFER! JOFFST+JJ.NFILE)  -  H(JJ) 

20  CONTINUE 

IRITER  •  MOD(NUMBERCNFILE) .NTBUF) 

C 

C  WRITE  BUFFER  SAMPLES  INTO  FILE  WHEN  BUFFER  IS  FULL 
C 

IF (IRITER  .EQ.  0)THEN 

WRITE(IUNIT)  NDATA1 . (RDATA1 (NN) ,NN»1 .NDATA1) 

WRITE (IUNIT)  NRLWRD. (BUFFER (KK.NFILE) ,KK»1 .NCXWRD) 
END  IF 
C 

C  WRITE  BUFFER  SAMPLES  INTO  FILE  WHEN  AT  END  OF  REALIZATION 
C 

IF(NUMBER(NFILE)  .EQ.  NTIMES  .AND.  IRITER  .NE.  0)THEN 
LFTCX  -  MOD (NTIMES, NTBUF) *NDELAT 
LFTRL  -  2*LFTCX 

WRITE(IUNIT)  NDATA1 . (RDATAl(NN) .NN-1 .NDATAl) 
WRITE(IUNIT)  LFTRL. (BUFFER (KK.NFILE) ,KK-1 .LFTCX) 
END  IF 
RETURN 
END 
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Table  A.9.  Listing  of  ACIRF  subroutine  READER. 


SUBROUTINE  READER(RESET , NDELAY , NTIMES . NFILE , IUNIT.H) 

C 

C*********************  VERSION  1.1  ***  08  JUL  1888  **************************** 

C 

C  THIS  SUBROUTINE  READS  A  COMPLEX  ARRAY  FROM  A  FILE  ON  UNIT-IUNIT 
C  ONE  DELAY  ARRAY  IS  READ  FOR  EACH  CALL  WITH  RESET  -  .FALSE. 

C  THE  INITIAL  HEADER  RECORDS  ARE  READ  WITH  RESET  -  .TRUE. 

C 

C  INPUTS  FROM  ARGUMENT  LIST: 

C 

C  RESET  -  LOGICAL  FLAG 

C  NDELAY  -  NUMBER  OF  DELAY  SAMPLES  IN  IMPULSE  RESPONSE  FUNCTION 
C  -  DIMENSION  OF  H  ARRAY 

C  NTIMES  *  NUMBER  OF  TIME  SAMPLES  IN  IMPULSE  RESPONSE  FUNCTION 
C  NFILE  -  FILE  NUMBER  (ONE  PER  ANTENNA) 

C  I UNIT  •  UNIT  OF  FILE 
C 

C  OUTPUTS  TO  ARGUEMENT  LIST: 

C 

C  H  -  IMPULSE  RESPONSE  FUNCTION  ARRAY  (DIMENSION  NDELAY) 

C 

C  COMMON  BLOCK  INPUTS: 

C 

C  RDATA1  *  HEADER  RECORD 
C 

C  PARAMETERS: 

C 

C  MAXBUF  «  MAXIMUM  BUFFER  SIZE  (RECORD  SIZE  IN  FILE  -  MAXBUF  +  1) 

C  MTENNA  -  MAXIMUM  NUMBER  OF  ANTENNAS  (EQUAL  TO  THE  MAX  NUMBER  OF  FILES) 

C  NDATA1  «  SIZE  OF  HEADER  RECORD  1 

C  NDATA2  -  SIZE  OF  HEADER  RECORD  2 

C 

O**************************************************************************** 

s 

PARAMETER  (MAXBUF-4086 . MTENNA-8 , NDATA1 -25 , NDATA2-32 ) 

COMMON/HEAD ER/RDATA1  (NDAl'Al) .  RDATA2  (MTENNA .  NDATA2) 

COMPLEX  H( NDELAY) . BUFFER (MAXBUF/2, MTENNA) 

DIMENSION  HEADER(NDATAl) .NUMBER(MTENNA) 

LOGICAL  RESET 

: 

:**♦**♦  IF  RESET  -  .TRUE.  READ  HEADER  RECORDS  AT  BEGINNING  OF  FILE 

: 

IF (RESET) THEN 
REWIND  I UN IT 
READ(IUNIT) 

READ(IUNIT)  MDATA1 , (HEADER(II) , II-l ,MIN0(MDATA1 .NDATAl) ) 

READ(IUNIT) 

MDATA1  -  MIN0(MDATA1, NDATAl) 


A- 14 


Table  A.9.  Listing  of  ACIRF  subroutine  READER  (Continued). 


DO  10  NN-1.MDATA1 
C 

C  FATAL  ERROR  IF  HEADER  RECORD  DOES  NOT  AGREE  WITH  RDATA1 
C  IF(HEADER(NN)  .HE.  RDATA1 (NN) )THEN 

WRITE! * .4000) NN , HEADER (NN) . UK , RDATA1 (KN) 

4000  FORMAT (21. 'FATAL  ERROR  REREADING  HEADER  RECORD  \/2X. 

1  'WORD  1 ,12.*  OF  HEADER  RECORD  IS  ’.1PE12.6/2X. 

2  ’WORD  1 .12, ‘  OF  HEADER  RECORD  SHOULD  BE  ’.E12.6) 

STOP 

END  IF 

10  CONTINUE 

NTBUF  -  !MAXBUF/2)/NDELAY 
MAXCX  -  NTBUF*NDELAY 
DO  20  NN-l.MTENNA 
NUMBER(NN)  -  0 
20  CONTINUE 

RESET  -  .FALSE. 

RETURN 
END  IF 

NUMBER(NFILE)  -  NUMBER(NFILE)  ♦  1 
C 

C  COUNT  NUMBER  OF  CALLS  TO  PROGRAM 

C  FATAL  ERROR  IF  ATTEMPT  TO  READ  BETOND  END  of  FILE 

C 

IF (NUMBER(NFILE)  .GT.  NTIMES)THEN 

WRITE!* ,4001) NUMBER (NFILE) , NTIMES 

4001  F0RMAT(2X,  'ATTEMPT  TO  READ  BETOND  END  of  FILE\/2X, 

1  'TIME  INDEX  •  ’ ,  I10/2X ,  ’MAX  INDEX  -  M10) 

STOP 
END  IF 
C 

C  READ  DATA  INTO  A  BUFFER 
C 

IF !MOD!NUMBER!NFILE) , NTBUF)  .EQ.  1)THEN 

READ(IUNIT)  MDATA1 , (HEADER(NN) .NN-1 ,MIN0!MDATA1 .NDATA1) ) 
READ(IUNIT)  NRLWRD , 

1  !BUFFER!KK,NFILE) ,KK«1 ,MIN0!NRLWRD/2,MAXBUF/2)) 

MDATA1  -  MINO(MDATAl.NDATAl) 

DO  30  NN-1.MDATA1 
C 

C  FATAL  ERROR  IF  HEADER  RECORD  DOES  NOT  AGREE  WITH  RDATA1 
C 
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Table  A.9.  Listing  of  ACIRF  subroutine  READER  (Concluded) 


IF(HEADER(NN)  .HE.  RDATA1(NN))THEN 

WRITE (* ,4000)NN,HEADER(NN) .  NN.RDATAl(NN) 
STOP 
END  IF 

30  CONTINUE 
END  IF 
C 

C  READ  BUFFERED  DATA  INTO  ARRAT 

3 

JOFFST  «  M0D( (NUMB ER(NFILE)-1)*ND ELAT. MAXCX) 

DO  40  JJ-l.NDELAT 

H(JJ)  -  BUFFER( J0FFST+ JJ , NFILE) 

40  CONTINUE 
RETURN 
END 
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The  subroutine  READER  reverses  this  process.  It  first  reads  a  record  into  a 
buffer,  and  then  outputs  a  single  impulse  response  function  array  from  the  buffer  each 
time  it  is  called.  After  (MAXBUF/2)/NDELAY  calls  to  READER,  another  record  is 
read  from  the  file  into  the  buffer,  and  so  on. 
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APPENDIX  B 


ACCURACY  OF  ANGULAR  INTEGRATION  TECHNIQUES 


B.l  INTRODUCTION 

The  mean  signal  energy  in  a  Kx  —  Kv  grid  cell  at  the  output  of  an  antenna 
is 


r(*+l/2)AK, 


E(k,l)  =  f 

J(k-l/2)AK, 


dK x 

2  7T 


/ 

Ju- 


(i+l/2)AA-y  dKv 


{l-l/2)AK,  27 r 


G{Kx-K0x,kv-K0v)Sk{Kx,Kv)  (B.l) 


This  equation  is  quite  general,  but  it  requires  that  the  energy  in  each  angular  bin  be 
calculated  and  stored  for  each  antenna.  This  latter  requirement  results  in  unaccept- 
ablely  large  arrays.  However,  Equation  B.l  can  be  approximated  by  assuming  that 
the  antenna  beam  profile  varies  slowly  over  the  grid  cells  so  G(K  —  K0)  may  be  pulled 
out  of  the  integral.  There  are  several  ways  that  this  can  be  done.  The  purpose  of  this 
appendix  is  to  calculate  the  accuracy  of  a  few  approximations  to  Equation  B.l. 

In  order  to  limit  the  scope  of  this  calculation,  isotropic  scattering  and  a 
uniformly  weighted  circular  antenna  will  be  assumed.  Without  loss  of  generality,  it 
can  then  be  assumed  that  the  antenna  is  pointed  away  from  the  line  of  sight  in  the  x 
direction,  or  equivalently,  that  the  pointing  azimuth  is  zero.  For  this  case,  the  angular 
part  of  the  GPSD  is 


Sk(Kx,K,}  =  tt£q  exp 


(ki + k ;k 

4 


(B.2) 


and  the  antenna  beam  profile  is 


G(KX,KV)  =  exp 


2  (g; + kik 

4 


(B.3) 


where  7  is  proportional  to  the  square  of  the  ratio  of  the  antenna  diameter  D  to  the 
decorrelation  distance: 
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4/n(2)  D 2 


(1.02899 tt)2  % 


2  * 


(B.4) 


For  this  isotropic  scattering  and  antenna  case,  the  power  at  the  output  of 
the  antenna  reduces  to 


P  1  i-xn  4/n(2)(^M>)J 
P a  =  — —  exp - — - 

1  +  Tf  1  +  7 


(B.5) 


where  rp  is  the  pointing  angle  and  80  is  the  beamwidth  of  the  antenna.  The  antenna 
filtered  decorrelation  distance,  which  is  the  same  in  both  the  x  and  y  directions  is 


f-A  —  \A  +  7^0 


(B.6) 


The  Kx  angular  grid  size  is  then 


A  Kx  = 


Nxtj\x  NxyJ  1  4-  7  A) 


(B.7) 


A  similar  expression  holds  for  A Kv. 

The  exact  expression  for  the  received  power  will  be  compared  to  the  total 
power  in  the  grid, 


Nz/ 2-1  JVy/2-1 

Pa=  £  £  E(k,l)  , 

k-~Nx!  2  l=-N,/ 2 


(B.8) 


to  compute  an  error  in  the  total  power 


PE  = 


\Pa~Pg\ 


(B.9) 


for  each  of  the  algorithms  used  to  evaluate  Equation  B.l. 
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B.2  ALGORITHMS 


The  first  approximation  to  the  exact  result  is  just  Equation  B.l.  This  ex¬ 
pression  will  result  in  some  error  in  the  total  power  because  of  round-off  errors  in  the 
summation  of  the  contributions  from  each  grid  cell  and  because  of  the  finite  size  of  the 
angular  grid.  Additional  error  results  when  the  beam  is  pointed  near  the  edge  of  the 
grid  so  part  of  the  beam  is  pointed  out  of  the  grid.  The  magnitude  of  this  error  will  be 
apparent  in  the  results  of  this  analysis.  With  the  assumptions  of  isotropic  scattering 
and  an  isotropic  Gaussian  beam  profile,  the  integrals  indicated  in  Equation  B.l  can 
be  obtained  in  closed  form  with  the  result 


A  (M  = 


%  1erfc 


erfc 


(k  -  l/2)AKzt<>y/TTl  1K0l0 


2  2y/lTl\ 

(k  +  1/2)  A  Kxto\/l  -f-  7 


2  vrn 


(B.10) 


erfc 


(/  -  1/2)A  Kytoy/l+^ 


—  erfc 


(/  +  l/2)AKv£o\/rT7 


where  erfc(-)  is  the  complementary  error  function  and  where 


K0to  =  (1.028997r)~^  .  (B.ll) 

"0 

The  second  approximation,  and  the  one  that  is  used  in  ACIRF,  is  obtained 
by  assuming  that  the  antenna  beam  profile  is  constant  over  an  angular  grid  cell  and 
can  therefore  be  pulled  out  of  the  integral.  The  result  is  the  product  of  the  antenna 
beam  profile  times  a  term  that  is  equal  to  the  incident  energy  in  a  grid  cell: 


EjM  =  G(K,  -  K0x,K,  -  K0.)  E,(k,l) 


(B.12) 


where  the  incident  energy  is 


E , 


(M)  =  f 


(k+  1/2)AK  j  dK  f{l+l/2)AK,  dKv 


l(k-l/2)&Kx  27T  J{l-1/2)AK ,  27T 

The  indicated  integrals  can  again  be  expressed  in  terms  of  error  functions: 


Sk(K.,Kv) 


(B.13) 
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E,(k,l) 


erfc 

(k  -  1/2)AKxIq 

—  erfc 

f(k  +  l/2)AifAl 

2 

[  2 

J 

erfc 

(/  -  l/2)AKyla 

2 

—  erfc 

(/  +  l/2)AKytc 
2 

11 

x 


(B.14) 


The  advantage  of  this  approximation  is  that  the  error  function  terms  depend  only  on 
the  environment  so  they  can  be  done  once  and  used  for  all  antennas  thereby  reducing 
the  required  processing  time  and  array  sizes. 

A  third  approximation  to  Equation  B.l  is  similar  to  the  previous  approxi¬ 
mation.  Rather  than  using  the  antenna  gain  at  the  center  of  each  angular  grid  cell, 
the  gain  averaged  over  the  grid  cell  is  used: 


E3(k,l)  =  Ej(k,l) 


(*+l/2)AJf,  r(i+l/2)AK, 


AKxAKy  J(k-l/2)AK: 


r 

j(k- 


iK‘  L. ... . .. '  <<*»  <3(Jf.  -  K»„  K ,  -  Ko.) 


/(t-l/2)AK, 


(B.15) 


Writing  the  indicated  integrals  in  terms  of  error  functions  gives  the  result: 


E3(k,l)  =  Et(k,l) 


7 r 


AKxAKy^ 


erfc 


erfc 


y/i\{k-l/2)AKx-K0]lo 

2 


y/i[{k  +  1/2)  A Kx  -  Ko 
2 


V7(/  -  l/2)AJfv£o 
2 


—  erfc 


-  erfc 

■Ji(i  + 1 [MKA 


1)  ■ 


I- 
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Finally,  a  simple  way  to  eliminate  the  problem  created  by  Equation  B.l  is 
to  do  away  with  the  integral.  This  zeroth  order  approximation  is 


A  k  A  K 

E4{k,1)  =  *  y  G{kAKx  -  K0tJAKv  -  K0y)SK{kAKx,lAKy)  .  (B.17) 
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B.3  RESULTS 


The  relative  error  (Equation  B.9)  of  the  four  algorithms  is  calculated  for  a 
range  of  the  ratio  of  the  decorrelation  distance  to  the  antenna  diameter.  The  scattering 
loss  of  the  antenna  for  this  range  of  Lo/D  is  shown  in  Figure  2.5  for  pointing  angles  of  0, 
0o/2,  and  0O.  Figures  B.l,  B.2  and  B.3  show  the  relative  errors  of  the  four  algorithms 
for  the  same  set  of  pointing  angles.  For  these  calculations,  a  32x32  angular  grid  is 
used.  This  is  the  smallest  angular  grid  size  used  in  ACIRF. 

When  the  pointing  angle  is  zero,  algorithm  1  (solid  line)  has  a  relatively 
constant  error  of  about  2  x  10  -5.  This  error  is  attributed  to  round-off  errors  that 
result  in  performing  the  sum  of  E(k,l)  over  all  grid  cells.  Both  algorithms  2  (dashed 
line)  and  3  (dotted  line)  for  this  case  have  peak  errors  of  about  10-3  which  occur  when 
the  angular  spread  of  the  incident  energy  is  about  equal  to  the  antenna  beamwidth  (i.e. 
£o  is  approximately  equal  to  D ).  Algorithm  4  (dash-dot-dash  line)  has  a  maximum 
relative  error  for  large  values  of  lo/D  of  about  2.5,  which  is  unacceptably  large. 

As  the  pointing  angle  increases,  the  maximum  error  of  algorithm  1  increases 
due  to  the  fact  that  part  of  the  beam  is  pointing  out  of  the  angular  grid.  However, 
even  when  the  pointing  angle  is  equal  to  a  beamwidth,  this  error  is  less  than  2  percent, 
and  it  is  less  than  0.1  percent  when  the  pointing  angle  is  one-half  of  a  beamwidth. 
Algorithms  2  and  3  generally  have  errors  which  are  close  to  that  of  algorithm  1.  The 
maximum  error  of  these  algorithms  is  less  than  2.5  xlO-3  when  the  pointing  angle 
is  0o/2.  When  the  pointing  angle  is  equal  to  a  beamwidth,  the  sign  of  the  errors  of 
algorithms  2  and  3  change  over  the  range  of  lo/D  resulting  in  errors  that  are  smaller 
in  magnitude  than  that  of  algorithm  1. 

Algorithm  2  generally  has  a  slightly  smaller  error  than  algorithm  3  over  the 
ranges  of  Iq/D  and  pointing  angles  considered  in  this  analysis.  Algorithm  2  is  also 
easier  to  implement  than  is  algorithm  3.  Based  on  these  results,  the  second  algorithm 
is  used  in  ACIRF. 
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Figure  B.2.  Relative  error  tor  pointing  angle  of  Qq/2. 
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lo/D 


Figure  B.3.  Relative  error  for  pointing  angle  of  9{ 
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ATTN:  IGARFUNKEL 
ATTN:  J  KLUCK 
ATTN:  MROLENZ 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
ATTN:  SECURITY 

ATLANTIC  RESEARCH  SERVICES  CORP 
ATTN:  R  MCMILLAN 

ATMOSPHERIC  AND  ENVIRONMENTAL  RESEARCH  INC 
ATTN:  MKO 

AUSTIN  RESEARCH  ASSOCIATES 
ATTN:  J  THOMPSON 

AUTOMETRIC  INCORPORATED 
ATTN:  C LUCAS 

BDM  INTERNATIONAL  INC 
ATTN:  L JACOBS 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  C  PRETTIE 
ATTN:  J  WORKMAN 
ATTN-  S BRECHT 

BOEING  CO 

ATTN:  GHALL 

CALIFORNIA  RESEARCH  &  TECHNOLOGY,  INC 
ATTN:  M  ROSENBLATT 

CHARLES  STARK  DRAPER  LAB,  INC 
ATTN:  ATETEWSKI 

COMMUNICATIONS  SATELLITE  CORP 
ATTN:  G  HYDE 

CORNELL  UNIVERSITY 

ATTN:  D  FARLEY  JR 
ATTN:  M  KELLY 

ELECTROSPACE  SYSTEMS,  INC 
ATTN:  P  PHILLIPS 


EOS  TECHNOLOGIES,  INC 
ATTN:  B  GABBARD 
ATTN:  W  LELEVIER 

GENERAL  ELECTRIC  CO 
ATTN:  REDSALL 

GENERAL  RESEARCH  CORP  INC 
ATTN:  J  EOLL 

GEO  CENTERS,  INC 

ATTN:  E  MARRAM 

GRUMMAN  AEROSPACE  CORP 
ATTN:  J  DIGLIO 

GTE  GOVERNMENT  SYSTEMS  CORPORATION 
ATTN:  W  I  THOMPSON,  III 

HARRIS  CORPORATION 
ATTN:  EKNICK 

HSS,  INC 

ATTN:  D HANSEN 

INFORMATION  SCIENCE,  INC 
ATTN:  W  DUDZIAK 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN:  E  BAUER 
ATTN:  H  WOLFHARD 

JS  LEE  ASSOCIATES  INC 
ATTN:  DR  J  LEE 

JAYCOR 

ATTN:  J  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
ATTN:  CMENG 
ATTN:  J  D  PHILLIPS 
ATTN:  R  STOKES 
ATTN:  T  EVANS 

KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 
ATTN:  G  DITTBERNER 

KAMAN  SCIENCES  CORPORATION 
ATTN:  B  GAMBILL 
ATTN:  DASIAC 
ATTN:  R RUTHERFORD 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  HENLEY 
ATTN:  JKUMER 
ATTN:  R  SEARS 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  DKREJCI 

LTV  AEROSPACE  &  DEFENSE  COMPANY 
2  CYS  ATTN:  LIBRARY 
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M  I  T  LINCOLN  LAB 

ATTN:  D  TOWLE  L-230 
ATTN:  I  KUPIEC  L-100 
ATTN:  MLEEG91 

MARTIN  MARIETTA  DENVER  AEROSPACE 
ATTN:  H  VON  STRUVE  III 
ATTN:  50018  J  BENNETT 

MAXIM  TECHNOLOGIES,  INC 
ATTN:  BRIDGEWAY 
ATTN:  J  SO 

MCDONNELL  DOUGLAS  CORP 
ATTN:  TCRANOR 

MCDONNELL  DOUGLAS  CORPORATION 
ATTN:  J  GROSSMAN 
ATTN:  RHALPRIN 

METATECH  CORPORATION 
ATTN:  R  SCHAEFER 
ATTN:  WRADASKY 

METEOR  COMMUNICATIONS  CORP 
ATTN:  R  LEADER 

MISSION  RESEARCH  CORP 
ATTN:  J  KENNEALY 
ATTN:  R  ARMSTRONG 
ATTN:  R  LARKIN 
ATTN:  S  BELANGER 
ATTN:  W  WHITE 

MISSION  RESEARCH  CORP 
ATTN:  BR  MILNER 
ATTN:  CLONGMIRE 
ATTN:  D  ARCHER 
ATTN:  DKNEPP 
ATTN:  D  LANDMAN 
ATTN:  FFAJEN 
ATTN:  F  GUIGLIANO 
ATTN:  G  MCCARTOR 
ATTN:  KCOSNER 
ATTN:  M  FIRESTONE 
ATTN:  RBIGONI 
ATTN:  R  BOGUSCH 
2  CYS  ATTN:  R  DANA 

ATTN:  R  HENDRICK 
ATTN:  RKILB 
ATTN:  S  GUTSCHE 
ATTN:  TECH  INFO  CENTER 
ATTN:  TECH  LIBRARY 

MITRE  CORPORATION 

ATTN:  D  RAMPTON,  PH.D 
ATTN:  MRDRESP 

MITRE  CORPORATION 

ATTN:  J  WHEELER 
ATTN:  MHORROCKS 
ATTN:  RCPESCI 
ATTN:  W  FOSTER 

NORTHWEST  RESEARCH  ASSOC,  INC 
ATTN:  E  FREMOUW 


PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  E  FIELD  JR 
ATTN:  F  THOMAS 
ATTN:  HBRODE 

PHOTOMETRICS,  INC 

ATTN:  ILKOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATTN:  D  BURWELL 
ATTN:  0  LEWIS 

PHYSICAL  RESEARCH  INC 
ATTN:  W.SHIH 

PHYSICAL  RESEARCH  INC 
ATTN:  HFITZ 
ATTN:  PLUNN 

PHYSICAL  RESEARCH,  INC 
ATTN:  R  DELIBERIS 
ATTN:  T  STEPHENS 

PHYSICAL  RESEARCH,  INC 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  W  SCHLUETER 

R  &  D  ASSOCIATES 

ATTN:  C  GREIFINGER 
ATTN:  F  GILMORE 
ATTN:  GHOYT 
ATTN:  M  GANTSWEG 

R&D  ASSOCIATES 

ATTN:  J  WALTON 

RAND  CORP 

ATTN:  C  CRAIN 
ATTN:  E  BEDROZIAN 

RAND  CORP 

ATTN:  B BENNETT 

RJO  ENTERPRISES/POET  FAC 
ATTN:  A  ALEXANDER 
ATTN:  W BURNS 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  T  CARNEY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  S  ROSENCWEIG 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  C  SMITH 
ATTN:  D  HAMLIN 
ATTN:  D SACHS 
ATTN:  ESTRAKER 
ATTN:  LLINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  COCKAYNE 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DTELAGE 
ATTN:  M  CROSS 
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SRI  INTERNATIONAL 

ATTN:  D  MCDANIEL 
ATTN:  J  DEPP 
ATTN:  R LEONARD 
ATTN:  W  CHESNUT 
ATTN:  WJAYE 

STEWART  RADIANCE  LABORATORY 
ATTN:  RHUPPI 

TELECOMMUNICATION  SCIENCE  ASSOCIATES 
ATTN:  R  BUCKNER 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  J  WOLFSBERGER,  JR 
ATTN:  N  PASSINO 

TOYON  RESEARCH  CORP 
ATTN:  J  ISE 

TRW  INC 

ATTN:  DRDGRYBOS 
ATTN:  R  PLCBUCH 
ATTN:  R5/220  H  CULVER 


TRW  SPACE  &  DEFENSE  SYSTEMS 
ATTN:  DM  LAYTON 

USER  SYSTEMS,  INC 

ATTN:  S  W  MCCANDLESS,  JR 

UTAH  STATE  UNIVERSITY 
ATTN:  K BAKER 

ATTN:  L  JENSEN,  ELEC  ENG  DEPT 

VISIDYNE,  INC 

ATTN:  J CARPENTER 

FOREIGN 

FOA  2 

ATTN:  B  SJOHOLM 

FOA  3 

ATTN:  T  KARLSSON 
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